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Abstract 

We study the correction-to-scaling exponents for the two-dimensional self- 
avoiding walk, using a combination of series-extrapolation and Monte Carlo meth- 
ods. We enumerate all self-avoiding walks up to 59 steps on the square lattice, and 
up to 40 steps on the triangular lattice, measuring the mean-square end-to-end 
distance, the mean-square radius of gyration and the mean-square distance of a 
monomer from the endpoints. The complete endpoint distribution is also calcu- 
lated for self-avoiding walks up to 32 steps (square) and up to 22 steps (triangular). 
We also generate self-avoiding walks on the square lattice by Monte Carlo, using 
the pivot algorithm, obtaining the mean-square radii to ~ 0.01% accuracy up to 
N = 4000. We give compelling evidence that the first non-analytic correction 
term for two-dimensional self-avoiding walks is Ai = 3/2. We compute several 
moments of the endpoint distribution function, finding good agreement with the 
field-theoretic predictions. Finally, we study a particular invariant ratio that can 
be shown, by conformal-field-theory arguments, to vanish asymptotically, and we 
find the cancellation of the leading analytic correction. 



KEY WORDS: Self- avoiding walk; polymer; exact enumeration; series expansion; 
Monte Carlo; pivot algorithm; corrections to scaling; critical exponents; conformal in- 
variance. 



1 Introduction 



The study of the universal properties of self- avoiding walks (SAWs) in the long-chain 
limit has been a central problem in both statistical mechanics and polymer physics 
for more than three decades [1-3]. In an iV-step chain, the mean value of any global 
observable O typically has an asymptotic expansion as N — > oo of the form 



(0) N = AN p < 
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where the leading exponent po and the correction-to-scaling exponents Ai < A 2 < ■ ■ ■ 
are universal, i.e. they depend on the spatial dimension d but not on the microscopic 
details of the interactions (provided that the interactions are short-ranged and primar- 
ily repulsive). This universality justifies the intense efforts that have been devoted to 
determining these universal exponents, using a variety of analytical and numerical ap- 
proaches. 

In this paper we will address the problem of determining the leading non-analytic 
correction-to-scaling exponent Ai for the two-dimensional self-avoiding walk and for 
the closely related problem of self-avoiding polygons (SAPs). At least two different 
theoretical predictions have been made for the purportedly exact value of this exponent: 
Ai = 3/2 based on Coulomb-gas arguments [4,5], and A x = 11/16 based on conformal- 
invariance methods [6]. In addition, several numerical methods have been employed to 
estimate Ai, notably exact enumeration and extrapolation (series analysis) [7-19] and 
Monte Carlo [18,20-23]. The estimates of Ai resulting from these numerical works are, 
for the most part, wildly contradictory: even when one compares estimates produced by 
a single method, such as series analysis, they range from ps 0.5 [18] to ~ 0.65 [8,10,13,14] 
to ~ 0.85 [15] to ~ 1 [7,14] to 1.5 [16,17,19]. Similar variation can be found in estimates 
of Ai obtained from Monte Carlo studies, ranging from ps 0.5 [18] to ps 0.6 [23] to 
ps 0.84 [21] to ps 1.1 [23] to ^ 1.2 [20]. 

Other models in the same universality class have also been considered, yielding results 
in contrast with those for the SAW. For instance, for lattice trails (connected paths where 
the self-avoidance constraint is applied only to bonds, not vertices) it was shown by a 
transfer-matrix study [24] that the correction-to-scaling exponent is indeed ~ 11/16, 
confirming an earlier result based on series analysis [25]. This same transfer- matrix 
study also found Ai ~ 3/2 for SAWs. What remains to be understood is why the 
contribution with Ai = 11/16 seems to be absent for square- and triangular-lattice 
SAWs, yet present for trails. 

To further confuse the subject, we should mention the recent results of Jensen [26] 
for osculating SAPs on the square lattice: these are a superset of SAPs in which bonds 
may touch at a vertex but not cross. In a sense they interpolate between SAPs, in which 
intersections are strictly forbidden, and closed trails, in which crossing at a vertex is 
allowed. A careful analysis of the corresponding series [26] shows very convincingly that 
Ai = 3/2 and finds no correction corresponding to Ai = 11/16. 
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To add to our list of unexplained phenomena, we remark on a recent Monte Carlo 
study of SAWs on the Manhattan lattice [27], where it was found that the critical 
exponent 7 has the same value 43/32 as for SAWs on regular lattices provided that the 
value Ai = 11/16 is used in the analysis. It is quite unclear why a different value of Ai 
should arise for the Manhattan lattice than is found for the square lattice. 

Returning now to the simplest case of SAWs and SAPs on regular lattices, the gross 
disparities among the extant estimates of the correction-to-scaling exponent might lead 
one to suspect that different methods are computing different quantities. For example, 
it might be that some methods are measuring (or predicting) the leading correction ex- 
ponent Ai, while others are measuring (or predicting) a subleading correction exponent 
A 2 or A 3 , and still others are measuring some sort of "effective" exponent A e ff that rep- 
resents phenomenologically the observed corrections to scaling in some specified interval 
of walk length N (and arising in reality from the sum of two or more correction-to-scaling 
terms) . 

Moreover, it may even be true that different observables produce different patterns of 
nonvanishing corrections to scaling. For instance, the Ai = 3/2 correction term appears 
to be present for SAWs but absent for SAPs. While this may at first sight be considered 
a violation of universality, we show below that it is not. 

Two recent analyses [17, 19] based on very long series for square-lattice SAWs and 
SAPs have, however, yielded a consistent and convincing picture of the corrections to 
scaling: the first non-analytic correction-to-scaling exponent is indeed just Ai = 3/2, as 
predicted by Nienhuis [4,5]; but there are also analytic corrections to scaling proportional 
to integer powers of 1/N, the first of which dominates asymptotically. More precisely, 
a careful numerical study based on a 51-term SAW series [17] found that the number of 
iV-step SAWs on the square lattice is given asymptotically by 

c N ~ /i JV A^ 11/32 [1.177043 + 0.5500/A^-0.140/A^ 3/2 -0.12/A^ 2 + ---] 

+ (-{if N~ 3/2 [-0.1899 + 0.175/iV - 1.51/N 2 + •••]. (1.2) 

It is likely that previous studies identified some sort of effective exponent that reflects 
a combination of the effects of the 1/N and 1/iV 3 / 2 correction-to-scaling terms (see 
Section 2.2 for further discussion of this point). Similarly, a careful numerical study 
based on a 90-term SAP series [19] found that the number of 2iV-step SAPs on the 
square lattice is given asymptotically by 

P2N ~ I2 2N N~ 5/2 [0.0994018 - 0.02751/N + 0.0255/N 2 + 0.12/N 3 + •••]; (1.3) 

note that here there is no N~ 3 ^ 2 (or N~ 5 ^ 2 ) correction term. Finally, a recent transfer- 
matrix analysis [24] of SAWs on the square lattice also found compelling numerical 
evidence in favour of the value Ai = 3/2, and against all values Ai significantly less 
than 3/2. 

One possible cause of some confusion is that, because the value of the leading critical 
exponent of the SAP generating function is 2 — a = 3/2, any correction-to-scaling term 
with A = half-integer "folds into" the analytic background term and is therefore unde- 
tectable! In other words, no corrections proportional to iV~ A appear in the coefficients 
Pn- In order to understand this point, let us recall that the critical exponent a is defined 
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by the leading asymptotic behavior p^ oc n N N a ~ 3 of the polygon counts, corresponding 
to a leading behavior P(x) = J2n>oPnx n ~ const x (1 — x/x c ) 2 ~ a as x | a; c = 
for the polygon generating function. If we now include both analytic and non-analytic 
corrections to scaling, the polygon generating function can be written generically as 

P(x) = J2p N x N ~ A(x) + B(x)(l-x/x c ) 2 ~ a [l + c(l-x/x c ) Al + ■■■} (1.4) 

Af>0 

with A(x) and B(x) analytic in the neighbourhood of the critical point x c = 1/ ji. Since 
a = 1/2 for two-dimensional SAPs [4,5,19], if A 1 = half-integer the above equation may 
be rewritten as 

P(x) = Y.P2nx N ~ A(x)+B(x)(l-x/x c ) 3 ' 2 [l + ---], (1.5) 

JV>0 

with the (1 — x/x c ) Al correction term absorbed into the analytic background term A(x). 
Therefore, if (1.4) holds, no correction of the form N~ Al is present. On the other hand, 
if the polygon counts were to exhibit a behaviour of the form p N oc fi N N a ~ 3 [l + . . . + 
a/N Al + . . .] with a = 1/2 and A x = 3/2 — and hence include a term oc ji N N~ A 
- then the generating function P(x) would exhibit, on top of the (1 — x/x c ) 3 ^ 2 leading 
behaviour, a non-analytic confluent term of the form (1 — x/x c ) 3 \og(l—x/x c ) in addition 
to the analytic term (1 — x/x c ) 3 . However, as discussed in some detail in [19], there is 
no evidence for a term of the form a/N 3 ^ 2 in the analysis of the SAP count series, and 
indeed there is considerable evidence for the absence of such a term. There is, however, 
abundant evidence of such a term in the radius-of-gyration series of SAPs [28]. 

In the present paper we make some further progress in supporting the assertion that 
Ai = 3/2 for SAWs on regular two-dimensional lattices (here square and triangular). 
First, we make a conventional analysis of corrections to scaling in the standard observ- 
ables (R 2 ), (R 2 ,) and (R^)', our contribution here is to present and use extended series 
expansions and a more efficient Monte Carlo algorithm. The results of this analysis 
are consistent with other recent work in supporting the conclusion that A 1 = 3/2. In 
the course of this analysis we point out that, for certain observables, pairs of correction 
terms of opposite sign can (and do) conspire to give an effective exponent that is smaller 
than both of the individual exponents; this explains the apparent exponents A < 1 
observed in some earlier work. Second — and this is perhaps our main contribution — we 
point out several observables in which a correction-to-scaling term becomes the leading 
term. These include: (a) the combination ^(R 2 ,) — 2(-R„) + ^(R 2 ), which arises in 
the conformal-invariance theory [29,30]; and (b) quantities related to the breaking of 
Euclidean invariance down to the lattice symmetry group, the simplest of which is (on 
the square lattice) the fourth-order moment (r 4 cos4#) = (x 4 — 6x 2 y 2 + y A ). Analysis 
of these quantities by Monte Carlo methods yields only a modest improvement over the 
analysis of conventional quantities — the trouble is that the new quantities exhibit a low 
"signal-to- noise ratio" — but the series analysis is quite precise. 

The plan of this paper is as follows. In Section 2 we define the quantities to be 
studied and collect some theoretical results that will be used or tested in the following 
sections. Section 3 reports the results of our series analysis: first, we analyze the SAW 
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counts (Section 3.3); then, we analyze the radius of gyration, the end-to-end distance 
and the average distance of a monomer from the endpoints, along with their invariant ra- 
tios (Section 3.4); finally, we analyze the higher-order rotationally-invariant moments of 
the endpoint distribution function (Section 3.5) and the corresponding non-rotationally- 
invariant moments (Section 3.6). For each of them, we determine the asymptotic be- 
haviour as N — > oo, focusing in particular on the correction-to-scaling exponent Ai 
and on the behaviour at the antiferromagnetic singularity (in the case of the square 
lattice). In Section 4 we report the analyses of our Monte Carlo data, confirming the 
absence of a correction-to-scaling exponent A 1 = 11/16. Finally, in Section 5 we draw 
our conclusions. 



2 Definitions and theoretical background 

In this section we review briefly the basic facts and conjectures about the SAW that will 
be used (or tested) in the remainder of the paper. 

2.1 Definitions and notation 

Let £ be some regular <i-dimensional lattice. Then an N-step self- avoiding walk (SAW) 
uj on C is a sequence of distinct points u>o, u>i, ujn in C such that each point is a nearest 
neighbour of its predecessor. We assume all walks to begin at the origin {ojq = 0) unless 
stated otherwise. 

First we define the quantities relating to the number (or "entropy") of SAWs. Let 
cn [resp. cjv(x)] be the number of iV-step SAWs on C starting at the origin and ending 
anywhere [resp. ending at x]. Then cn and Cjv(x) are believed to have the asymptotic 
behaviour 

c N ~ const x ^N 1 - 1 (2.1) 
c N (x) ~ const x fi N N a - 2 (x fixed ^ 0) (2.2) 

as iV — > oo; here \i is called the connective constant of the lattice, and 7 and a are 
critical exponents. The critical exponents are believed to be universal among lattices 
of a given dimension d. For rigorous results concerning the asymptotic behaviour of cat 
and cn(x), see [31-34]. 

Next we define several measures of the size of an iV-step SAW: 

• The squared end-to-end distance 

Rl = <*>%■ (2-3) 

• The squared radius of gyration 

Z-wn^r £}<*-»>)•■ (") 

• The mean-square distance of a monomer from the endpoints 
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R 



2(N 



N 

■yyEh 2 + (^-^v) 2 ] • (2-5) 



We then consider the mean values (R 2 )n, {R 2 )n and (Rf^N hi the probability distribu- 
tion that gives equal weight to each iV-step SAW. Very little has been proven rigorously 
about these mean values, but they are believed to have the leading asymptotic behaviour 

(R 2 e ) N , (R 2 g ) N , (R 2 J N ~ const x N 2 » (2.6) 

as N — > oo, where v is another (universal) critical exponent. Hyperscaling [35] predicts 
that 

dv = 2- a. (2.7) 

For SAWs in two dimensions, Coulomb-gas arguments [4, 5] as well as arguments 
based on stochastic Loewner evolution (SLE) [36] predict that u — 3/4, a — 1/2 and 7 = 
43/32. Prior numerical studies have confirmed these values to high precision [17,19,37]; 
in this paper we take them for granted. 

The amplitude ratios 

(Rg)N . , 

\Re)N 

R _ (Rmh f q\ 

Bn - im\ ( 2 - 9 > 

\Re)N 

are expected to approach universal values in the limit N — > 00, which we call A and 
B; one of our goals is to estimate these limiting amplitude ratios. Many other universal 
amplitude combinations (notably involving SAPs) are discussed in [38,39]. 
Of particular interest is the linear combination [29,30] 

(2.10) 




Vh, 

and the corresponding unnormalized quantity 



f N = F N (R 2 e ) N = [2 + (R 2 g ) N - 2(R 2 J N + l -(R 2 e ) N , (2.11) 

where y t — 1/v and = 1 +7/ (2u) are the thermal and magnetic renormalization-group 
eigenvalues, respectively, of the n- vector model at n — 0. In two dimensions — where 
y t = 4/3 and y h = 91/48, hence 2 + y t /y h = 246/91 — Cardy and Saleur [29] (as corrected 
by Caracciolo, Pelissetto and Sokal [30]) have predicted, using conformal field theory, 
that limjv^oo-^jv — 0. We shall henceforth refer to this relation as the CSCPS relation. 
This conclusion has been confirmed by previous high-precision Monte Carlo work [30] 
as well as by series extrapolations [40]. It is therefore of interest to examine the rate at 
which F N tends to zero, as this gives information on the correction-to-scaling terms. We 
will discuss this from a theoretical point of view near the end of Section 2.2, and from 
a numerical point of view in Sections 3.4.2 and 4.2. 
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It turns out that limjv->-oo Fn = holds not only for the ordinary square-lattice SAW, 
but also for SAWs with nearest-neighbour interactions, right up to (but not at) the theta 
point [41]. Moreover, the relation appears to hold at the theta point if 2 + y t /y h is given 
its theta-point value 23/8 instead of 246/91. This observation is used in [41] to locate 
the theta point more precisely. 

We shall also consider higher moments of the end-to-end distance. Limiting ourselves 
to two-dimensional lattices, let us write 

uj n = (x,y) = (r cos 9 , r sin 9) . (2-12) 

The Euclidean-invariant moments (r k ) N are of course expected to behave as 

(r k ) N ~ const x N ku (2.13) 

as N — > oo. One can consider the dimensionless ratios 

M2 k ,N = ^ , (2.14) 

V / N 

which approach finite limits for N — > oo; these limiting ratios M 2 fc i00 are universal quan- 
tities that characterize the end-to-end distribution function. Estimates of M 2 fc i00 have 
been obtained in [42] using field theory and the Laplace-deGennes transform method. 
It turns out [43, 44] that the 2-point function is very nearly equal to that of a free field, 
so that when the rescaled inverse propagator in momentum space 1 is written as 

oo 

DiqY 1 = l+g 2 + £M 2n , (2.15) 

n=2 

one has 1 > \b 2 \ > |6 3 | > \b 4 \ > . . . . One obtains [42] 

M2k,oo = n/ r( l +2 ^\ fc _ 1 [1 - h(k - 1) + R k \ k\ n (l + f ) , (2-16) 



T( 7 + 2A;z/)r(7)* 



3=0 



where Rk is a very small correction (unless k is very large) that involves the constants 
6 3 , 64, ... as well as higher powers bibj, bibjbk, .... Explicitly, 

k 1 
Rk = Y,(- 1 ) n+1 ( k + 1 ~ n ) b n + ^(k-2)(k-3)b 2 2 + ■■■ (2.17) 

n=3 

Note that R 2 = exactly. The universal nonperturbative constants fe 2 , fc 3 , . . . have 
been obtained from the analysis of exact-enumeration series on the square, triangular 
and hexagonal lattices [44]. Numerically, it is found [44] that b 2 is extremely small, 

1 D(q) is the Fourier transform of the two-point correlation function, rescaled so that the first two 
terms at small q are 1 — q 2 + 0(q 4 ). 
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b 2 = 0.00015(20), and that 63 is even smaller, |6 3 | < 3 x 10 5 . Using the estimate of 62 
in (2.16) and neglecting R k , we obtain for the lowest values of k: 

M 4 ,oc = 1.44574(28) (2.18) 

M 6j00 = 2.5876(10) (2.19) 

M 8j00 = 5.3805(32) (2.20) 

M 10iOO = 12.557(10). (2.21) 

A second class of interesting observables are moments that are invariant under the 
symmetry group of the lattice but not under the full Euclidean group: examples are 
the moments (r k cosn6) with n 7^ 0, where for the square (resp. triangular) lattice n 
must be a multiple of 4 (resp. 6). We expect these non-Euclidean-invariant moments to 
behave as 

(r k cos n9) ~ const x N kv ~ A ™ , (2.22) 

where A nr > is a new correction-to-scaling exponent [43] associated with the breaking 
of full rotation invariance down to the lattice rotation group: it thus depends on the 
lattice in question (e.g., square or triangular) and is in general different from the leading 
correction-to-scaling exponent Ai (which corresponds to a Euclidean-invariant irrelevant 
operator). 

For Gaussian models — and thus also for n-vector models (including the SAW case 
n = 0) in dimension d > 4 — we have A nr = 2v = 1 on any hypercubic lattice. For 
n-vector models in dimension d — 4 — e, this relation is modified at order e 2 [43]: 



A nr = v 



2 + + °< f3 > 



(2.23) 



In dimension d — 3, several alternative methods — field theory and exact-enumeration 
analysis — show that A nr is very close to 2v>, though not exactly equal [43]. In two 
dimensions on the square lattice, A nr = 2v exactly for the Ising model and for the n- 
vector model with n > 3 (in the latter case with logarithmic corrections) [43]. For the 
triangular lattice, similar arguments [43] predict A nr = 4z/. 2 For the Ising model, these 
predictions can be obtained using conformal field theory (see [45,46] for the classification 
of the subleading operators appearing in the Ising model); they can be checked explicitly 
[43] for at least one specific observable, using the analytic expression for the mass gap 
[44, 47, 48]. 3 It is therefore suggestive to conjecture that the same relations between 
A nr and v are valid for the SAW. This would predict A nr = 3/2 on the square lattice, 

2 For the hexagonal lattice, A nr = 4i/ for observables that break rotational invariance but are invariant 
under interchange of the two sublattices, while A nr = 3v for observables that distinguish the two 
sublattices [43]. 

3 For instance, consider on a square lattice the mass gap m(h) in the direction n defined as 



m(n) = - lim -i- log ( V G(x) ) , 

\x-h—r / 



were h = (cos 9, sin 9) is a unit vector and the summation runs over all x such that x ■ h — r. From the 
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and A nr = 3 on the triangular lattice. In Section 3.6 we will test (and confirm) this 
conjecture, by series analysis, for both square-lattice and triangular-lattice SAWs. 

2.2 Corrections to scaling 

Let us now make some general remarks concerning corrections to scaling. Clearly, 
(2.1)/(2.2)/(2.6) are only the leading term in a large- iV asymptotic expansion. Ac- 
cording to renormalization-group theory [49], the mean value of any global observable 
O behaves as N — > oo as 

(0) N = AN po 



Thus, in addition to "analytic" corrections to scaling of the form a,k/N k , there are 
"non-analytic" corrections to scaling of the form bk/N Al+k , Ck/N A2+k and so forth, 
as well as more complicated terms [not shown in (2.24)] which have the general form 
const/iV fclAl+fe2A2+ "' + ' where ki,k 2 ,-- - and I are non-negative integers. The leading 
exponent po and the correction-to-scaling exponents Ai < A 2 < . . . are universal; 
Po of course depends on the observable O in question, but the Aj do not. The various 
amplitudes (both leading and subleading) are all nonuniversal (and of course also depend 
on the observable 4 ). However, ratios of the corresponding amplitudes A, bo and Co (but 
not ak or the higher bk,Ck) for different observables are universal [50,51]. 

In fact, (2.24) is incomplete, as there are "mixing" terms arising from the fact that 
the temperature deviation from criticality is a smooth but nonlinear function of the 
nonlinear scaling fields g t and g^. This has the consequence [35, 52-55] that the sus- 
ceptibility (or SAW generating function), which has a leading singularity (x c — rc)~ 7 , 
also contains an additive term proportional to the energy, of order (x c — rr) 1- ". In the 
case of the two-dimensional Ising model, we have a — 0, and this term is responsi- 
ble for the logarithmic terms in the susceptibility, as was recently exhaustively studied 
in [56]. For the two-dimensional SAW, we have a = 1/2, and so one would expect a 
term A(x)(x c — x) 1 / 2 in the SAW generating function. To incorporate this term requires 
that the naively expected asymptotic form 

c N ~ fx N N ll / 32 [a + ai /N + a 2 /N 3 ' 2 + a 3 /N 2 + a A /N 5 ' 2 + • • •] (2.25) 

be modified to read 

c N ~ fM N N 11 / 32 [ao + a 1 /N + a 2 /N 3 / 2 + a 3 /N 2 + a 4 /N 5 / 2 + ---} 

+ fi N N~ 3/2 [S + a 1 /N + ■■■]. (2.26) 

exact solution [47], one can easily see that, for (3 — > /3 C , 

m(n) = mo(/3 c - /3)" 1 [l + (/? c - /3) 2 (a + fe o cos40) + • • •] 
with &o 7^ 0. This result shows explicitly that A nr = 2 = 2^. 

4 Sometimes a particular correction-to-scaling amplitude will vanish for some observables but not for 
others (e.g. for symmetry reasons). 



1 + ^ + ^ + 
N N 2 



c 



+ 



+ 



Cl 



_bo_ 

ATA, 



+ 



+ 



bi 
i\T A i +1 
c 2 



+ 



jV A i+ 2 



+ 



jyA 2 j\TA 2 +i j\r A2+2 



+ 



(2.24) 
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For loose-packed (i.e., bipartite) lattices, such as the square and simple-cubic lattices, 
there is an additional set of terms arising from the antiferromagnetic singularity, of the 
form 



l) N N q 



di d 2 e ei e 2 

N N 2 N A t F N A ? F+1 N A t F+2 



(2.27) 



where the exponent q of course depends on the observable. We know of no theoretical 
argument that predicts the value of the exponent A^ F . For the exponent q, in the 
closely related problem of the Ising-model susceptibility in two and three dimensions, 
Sykes [57] has given a configurational "counting theorem" that enables one to guess 
that the antiferromagnetic susceptibility behaves as the internal energy. This reasoning 
is discussed in greater detail in [58, 59]. 5 It follows that there should be a term in the 
susceptibility of the form D(x)(l + x/x c ) 1 ~ a [where D(x) is analytic in a neighbourhood 
of the antiferromagnetic critical point x = —x c ] and thus a term (—x c )~ N N a ~ 2 in the 
high-temperature-series coefficients. This result can be put on more solid ground [60] 
by noting that at the antiferromagnetic critical point the (unstaggered) magnetic field 
is an irrelevant variable, so that the leading contribution to the free energy is 

F(x,h) = ag t (x,h) 2 - a + F Teg (x,h), (2.28) 

where x is the inverse temperature and gt(x, h) is the nonlinear scaling field associated 
with the temperature at the antiferromagnetic critical point. Since 

g t (x,h) = (l + x/x c ) + a t h 2 + ... , (2.29) 

by performing the appropriate derivatives we obtain the result reported above (provided 
of course that a t ^ 0). This argument is very general and applies to any n- vector model; 
in particular, it applies for n = 0, i.e. to the SAW. Thus, for the SAW counts cn we 
expect a term {—fi) N N a ~ 2 , so that q = a — 2 for this observable [61]. 

The argument of Sykes can be generalized to higher moments of the two-point func- 
tion, i.e., J2 r \r\ 2k G(r). Also in this case one can identify two terms: one is proportional 
to the energy, while the other is conjectured to give an algebraically small (i.e., noncriti- 
cal) correction at the antiferromagnetic critical point. Such a conjecture was numerically 
verified in [62] for the three-dimensional Ising model. As for the susceptibility, this im- 
plies that asymptotically the moments have the form D(x)(l + :r/:r c ) 1-a , with D(x) 
analytic, for any k. Therefore, a term (— fi) N N a ~ 2 should be present in their high- 
temperature-series coefficients. Extending this conjecture to the n-vector model and in 
particular to the SAW (n — 0), we predict 

c N ~ ^N^iao + •••] + (-fi) N N a - 2 [d + ---} , (2.30) 
c N (r 2k ) N ~ ^N 2 ^- 1 [a' + •••] + (-fj,) N N a - 2 [d' + ---] . (2.31) 



5 The basic idea is that the susceptibility can be rewritten as the sum of two terms: one proportional 
to the energy, and a second one which can be argued (by series analysis) to give an algebraically small 
contribution near the antiferromagnetic critical point. 
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For k — 1, (2.31) gives the behaviour of cn(R 2 )n- It may seem natural to generalize 
the expression (2.31) to the other metric quantities, namely and R 2 g . Surprisingly 
(to us), our subsequent analysis (see Section 3.4.2) shows that, while the unnormalized 
second-moment series of the end-to-end distance series behaves precisely as expected in 
(2.31), the unnormalized series corresponding to both the radius of gyration and the 
mean monomer-endpoint distance behave a little differently. We find 

cn(RIJn ~ H N N^'\a' + •••] + {-n) N [d' Q + •••]. (2.32) 

That is to say, the antiferromagnetic exponent is different in the latter cases, namely 
instead of a — 2 = —3/2. Nevertheless, by taking the quotient of either (2.31) or (2.32) 
by (2.30), we obtain in all cases 

(R 2 ) N ~ 7V 2l H + ---] + (-l^TVK' + H (2-33) 

with q = 2u + a — 1— 7. [For the end-to-end distance, the dominant antiferromagnetic 
correction always comes only from (2.30); for the other two metric quantities, it comes 
from both (2.30) and (2.32).] For the end-to-end distance only, we have the additional 
relation 

Similarly, the rotationally-invariant higher moments (r 2k ) N are expected to behave 

as 

(r 2k ) N ~ N 2k »[a™ + ■ ■ ■} + (-l) N Ni*[d™ + ---} (2.35) 

with qt = 2kv + a — 1 — 7. The coefficients Oq" and d'{] also satisfy a relation analogous 
to (2.34). Our numerical analysis, described below, provides supporting evidence that 
the corresponding exponents are indeed qu = 3k/2 — 59/32 in two dimensions (see 
Sections 3.4.2 and 3.5). 

Finally, the non-analytic correction-to-scaling exponent A^ F was found numerically, 
in the case of the square-lattice SAW counts, to be 1 [17]. It would seem likely that this 
value should also hold for other properties, such as the metric quantities (R 2 )n- Our 
numerical studies, discussed below, are consistent with this conjecture — or, put another 
way, they are insufficiently sensitive to refute this obvious first guess. 

Let us now return to the question of the corrections to the CSCPS relation liniTv^oo Fn 
[cf. (2.10)]. Series analysis and Monte Carlo simulations (see Sections 3.4.2 and 4.2 
below) indicate that F N oc iV~ 3//2 , i.e. that the leading analytic correction cancels. Such 
a cancellation may seem surprising, but it can be understood by means of a standard 
renormalization-group argument. Consider the continuum 0{n) model with Hamiltonian 

H = H* + J d 2 r[tE(r) + hs 1 (r)] , (2.36) 

where H* is the fixed-point Hamiltonian, and E{r) and s l {r) are the energy and spin 
operators, respectively. The CSCPS relation [29, 30] is a consequence of the sum rule 

J d 2 r (6(0) cont e(r) cont ) = , (2.37) 
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where ©(r) cont is the trace of the continuum stress-energy tensor, and of course we must 
set n = to obtain SAWs. In order to translate this continuum relation into a lattice 
one, we must relate the continuum operator to its lattice counterpart. It is natural to 
assume that the trace of the lattice stress-energy tensor, 6(r) latt , whose explicit form is 
given in [29], behaves as 

0(r) latt = Z(t,h)Q(r) cont + ... , (2.38) 

where Z(t, h) is a smooth function of t and h, and the dots represent the contributions 
of the subleading operators. As a consequence of (2.38) we have 

J d 2 r (0(O) latt 6(r) latt ) = 0(t A \h^^ Al ) . (2.39) 

No corrections of order t appear in the previous relation. Eq. (2.39) therefore implies 
the absence of the analytic corrections in the CSCPS relation limjy^oo F N = 0. 

The observation that F N oc iV _3//2 implies a constraint on the subdominant ampli- 
tudes. More precisely, if we write 

(R 2 e ) N ~ a e N 3/2 + b e N 1/2 + c e + 0{1/ y/N) (2.40) 
(R 2 g ) N ~ a g N z / 2 + b g N 1 / 2 + c g + 0(1/ VN) (2.41) 
(R 2 Jn ~ a m N*' 2 + b m N 1 / 2 + c m + 0(1/ VN), (2.42) 

then the original CSCPS relation F N — > implies 

91a e = 364a m - 492a 9 , (2.43) 

while the absence of a 1/JV term in Fn means that the leading subdominant terms also 
satisfy an amplitude relationship analogous to (2.43), namely 

916 e = 3646 m - A92b g . (2.44) 

Note too from (2.11) that 

/ = lim f N = (2 + y T /y H )c e - 2c m + c g /2. (2.45) 

Let us conclude by discussing briefly the behaviour of "effective exponents" . Given 
a function f(N), let us define A eff (iV) by fitting locally to the Ansatz f (N) = a + b/N A : 
this gives 

Applying this to f(N) = a + ai/iV Al + a 2 /N A \ we obtain 

A (N) + ^A 2 N- A ^ 

Thus, if cii and a 2 have the same sign, the "effective" exponent A e ff(iV) lies between A x 
and A 2 for all N, and decreases monotonically to A x as — > oo. (This is the behaviour 
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one would expect intuitively for an "effective exponent".) By contrast, if a± and a 2 have 
opposite signs, then A eff (iV) starts above A 2 for small N, increases monotonically and 
reaches +00 at some finite value of N; it then jumps to —00, after which it continues 
to increase monotonically, tending asymptotically to A 1 from below as N — > 00. Thus, 
the qualitative behaviour of the effective exponents depends crucially on the sign (and 
magnitude) of ai/a 2 , which can vary from one observable to another. We shall see this 
phenomenon quite clearly in the two-dimensional SAW. 

3 Series analysis 

3.1 Summary of our data 

We have previously reported enumerations of square-lattice SAWs up to 29 steps for c^, 
(RI)n, (R 2 g )N and (R^n [63,64], and of triangular-lattice SAWs up to 22 steps for c N 
and (RI)n [63] and up to 19 steps for (R 2 )n and (R^n [64]. Analysis of these series can 
be found in [63,64]. We have also previously presented the square-lattice SAW counts cjv 
up to 51 steps [17] and the square-lattice polygon counts p 2 N up to 90 steps [19]. Analysis 
of these SAW series [17,19] provided good evidence that the non- analytic correction-to- 
scaling exponent is exactly Ai = 3/2 as predicted by Nienhuis [4,5], and that there is 
also the expected analytic term of leading order 1/N (as well as 1/N 2 , 1/iV 3 , . . .). For 
SAPs we found compelling evidence for purely analytic correction-to-scaling terms. We 
have thus far found no numerical evidence of a second non-analytic correction-to-scaling 
exponent A 2 , although it is reasonable to expect that one exists. 

In the present paper, we have extended the previous work by enumerating all SAWs 
on the square lattice up to 59 steps, and on the triangular lattice up to 40 steps, us- 
ing refinements of the finite-lattice method (FLM) due to Rogers (unpublished) and 
Jensen [65]. The results for cjv, (RI)n, (R 2 9 )n and (R^n are collected in Tables 1 and 
2. 

For square-lattice SAPs, the counts are now known up to 110 steps [66], and the 
radii of gyration up to 100 steps [66]. For triangular-lattice SAPs, the counts p N were 
previously known up to N = 35 [67]; in as- yet-unpublished work, one of us (ANR) has 
extended them up to iV = 40 [68], while even more recently another of us (I J) has 
extended the series to N = 60 steps [69]. 

As the FLM does not enable us to record the full end-point distribution, nor higher 
moments (at least not with the amount of memory available to us), we also programmed 
a conventional backtracking algorithm and recorded the full end-point distribution cat(x) 
for N < 32 (square lattice) and for N <22 (triangular lattice). As a result, we are able to 
study arbitrary moments. We refrain here from inundating the reader with the complete 
tables of Cat(x); they are available on www.ms.unimelb.edu.au/~iwan. However, we do 
list here most of those series that we subsequently analyse. For the square lattice, 
we give in Table 3 the rotationally invariant moments (r 4 )^, (r 6 ) N and (r 8 )^, and in 
Table 4 the corresponding non-rotationally-invariant moments (r 4 cos4#)jv, (r 6 cos4#)jv 
and (r 8 cos 46) n- For the triangular lattice, we give in Table 5 the rotationally invariant 
moments (r 4 )^, (r 6 )^ and (r 8 )^, and in Table 6 the corresponding non-rotationally- 
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invariant moments (r 6 cos6#)jv and (r 8 cos66>)at. 
3.2 Method of analysis 

In this subsection we explain in detail the method we used to analyse the data, the results 
of which are reported in subsequent subsections. For the triangular lattice, we expect 
the series coefficients of any generic quantity, such as the SAW generating function, to 
have an asymptotic expansion of the form 

^^{ao + t^j , (3.1) 

where \i is the connective constant and 7 is the critical exponent. Likewise, we expect 
the square-lattice series coefficients to have an asymptotic expansion of the form 

Z^ 7 " 1 («o + E jfc) + (~,) N N^ (b + ± ^p) , (3-2) 

where a is the critical exponent occurring in the polygon generating function. Similar 
expansions hold for metric quantities, and involve also the critical exponent v. Since the 
exact values 7 = 43/32, a — 1/2 and v = 3/4 are well established, we shall use them 
throughout this paper. 

Given the calculated terms of the series up to some order iV max , we proceed as follows: 
First we decide how many correction terms {a{} and {bi} we wish to include (i.e., we 
fix the numbers k and m); then we make some assumption for the values of //, A, and 
Af F ; finally, we fit the data to (3.1) or (3.2) by taking (k + m + 2)-tuples of successive 
values of N and solving for {a;} and {bi}. This can be done by solving a system of linear 
equations. 

By using (k + m + 2)-tuples at steadily larger values of N, many estimates for the 
{oj} and {h} are found. If the different estimates seem stable as N grows, we presume 
that they provide an acceptably accurate estimate of the actual asymptotic coefficients. 

A noteworthy feature of the method is that, if a blatantly-too-low correction-to- 
scaling exponent is given as input (for example, specifying A x = 1/2 for the two- 
dimensional SAW), the sequence of amplitude estimates for the term corresponding 
to that exponent will converge rapidly to zero, giving a very strong signal that the ex- 
ponent in question is absent. (Of course, if such a term were to occur with an amplitude 
several orders of magnitude smaller than the amplitudes of the other terms, one could 
be fooled into thinking such a term is absent. Our analysis assumes the absence of such 
pathologies.) 

Another point to bear in mind is that even if one knows the precise asymptotic 
form, with a limited number of series coefficients one can fit only to a small number of 
asymptotic terms (i.e., k and m cannot be taken too large). Beyond a certain number 
of terms in the asymptotic form, the quality of the fit visibly deteriorates. The more 
series coefficients are available, the more terms can be included in the Ansatz (provided 
that sufficient numerical precision is retained during the analysis). 
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A universally observed feature of the method is that the apparent accuracy of the 
amplitude estimates decreases rapidly as we move to higher-order terms in the asymp- 
totic expansions. That is to say, the apparent accuracy of the estimate of amplitude 
a i+i is significantly less than that of a; L . Moreover, adding further terms in the assumed 
asymptotic form (i.e., increasing k and m) improves convergence of the low-order am- 
plitudes di provided that k + m does not get too large, but after a certain point actually 
slows the convergence. In the case at hand, allowing more than 2-5 terms (these being 
the values of k + 1 and m + 1 separately) in the assumed asymptotic form led to a 
deteriorating (i.e., less stable) fit. 

As the series data at very small N are probably not reflective of asymptotic behaviour, 
and we have here the luxury of access to many terms (i.e., quite large N max ), the first 
19 — k — m terms of the series will not be used in any of our analyses here. 

Our analysis thus comprises two phases. In the first phase, we determine the cor- 
rect connective constant ji and the correct exponents Aj and Af F for the asymptotic 
expansion, as just described. In the second phase, we determine how many terms in 
the asymptotic expansion we can reliably use. We now describe our procedure for the 
second phase of the analysis. 

We begin by fitting for only one correction coefficient, a\. Then we add further 
asymptotic terms until the estimates obtained do not appear to be converging as iV — > 
oo to a value that is consistent with the previous estimates given by fits with one 
fewer asymptotic term. We define "consistent" by the requirement that estimates of all 
included asymptotic coefficients be well-converged and of the same sign and within a 
factor F = 2.4 of the previous estimates. More specifically, we invoke this requirement 
as follows: setting k = K gives estimates of a , . . . , a^; repeating the analysis with 
k — K + 1 yields estimates of a , . . . , a K+1 . We require that the coefficients a , . . . , clk 
from the two fits agree in sign and in magnitude within a factor of F; otherwise, we 
reject the fit with k — K + 1 and stop at k = K. Note that the choice of the value of F 
is somewhat arbitrary. Realistically, one can reasonably make any choice in the range 
1.5 < F < 3, the lower value being more conservative. We chose a value in this range 
that included most data sets with small values of k and m, and excluded those with 
higher values. 

Please note that the convergence (as N grows) of each fit is here judged by traditional 
intuitive (and thus somewhat subjective) methods. It would be an interesting project 
to find a precise definition of "well-converged" (or its synonym, "stable") that accords 
satisfactorily with our intuitive judgments and gives good results on test series; this 
would allow the series analysis to be converted into a precise algorithm. But we do not 
purport to carry out such a project here. 

For the triangular lattice, this procedure is thus relatively simple to implement. We 
compute fits initially with k — 0, incrementing A; by 1 until a non-stable or inconsistent 
estimate (as defined in the preceding paragraph) is found; we then revert to the previous 
group of stable and consistent estimates. The final entries (i.e., those corresponding to 
the maximum N) in the largest stable group are taken as our final estimates. 

For the square lattice, the procedure is more complicated, as it is not clear a priori 
whether terms involving Aj or Af F should be added to a given group. Empirically we 
have found that groups containing approximately equal numbers of Aj and Af- F terms, 
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or slightly more Aj terms, are more stable than estimates with significantly different 
values of k and m. Hence we begin by exploring groups with equal numbers of Aj and 
Af F terms, that is with k — m, adding one coefficient to each group at every stage. 
Next we try groups with one more Aj term than Af F terms, so that k — m + 1; and 
finally we try groups with two more Aj terms, so that k = m + 2. Again, the largest 
group that provides stable and consistent estimates is selected. As always, the given 
estimates are taken from the fits to the largest available value of N, which is N maX) since 
these should best reflect the asymptotic regime. 

The estimated error is calculated as the change between the estimate given by the 
longest series and the series ten terms shorter, multiplied by a factor reflective of the 
expected rate of convergence of the estimates. This latter factor is determined by as- 
suming that the error in the estimates is principally given by the first omitted Aj or 
Af F term. The difference in the exponents between the term in question and the first 
omitted term is then used to predict the value of the estimate on a fit to an infinite 
series. 

To illustrate these procedures, we show below the output from fitting the triangular- 
lattice series for Cn(R 2 )n with increasing numbers of correction-to-scaling terms. We 
make the Ansatz 

(N + l) 2 c N {R 2 g ) N /6 ~ f x N N l23 ^ 2 [a + a 1 /N + a 2 /N 3 / 2 + a 3 /N 2 + a i /N 5 / 2 + ---] (3.3) 
and obtain fits as follows: 



N 


a 


ai 




21 


0. 


01929438 


0. 


08276963 


22 


0. 


01932692 


0. 


08208623 


23 


0. 


01935556 


0. 


08145624 


24 


0. 


01938091 


0. 


08087327 


25 


0. 


01940346 


0. 


08033194 


26 


0. 


01942363 


0. 


07982768 


27 


0. 


01944175 


0. 


07935660 


28 


0. 


01945809 


0. 


07891535 


29 


0. 


01947289 


0. 


07850104 


30 


0. 


01948633 


0. 


07811114 


31 


0. 


01949859 


0. 


07774345 


32 


0. 


01950980 


0. 


07739603 


33 


0. 


01952007 


0. 


07706717 


34 


0. 


01952952 


0. 


07675536 


35 


0. 


01953823 


0. 


07645924 


36 


0. 


01954628 


0. 


07617761 


37 


0. 


01955373 


0. 


07590938 


38 


0. 


01956064 


0. 


07565359 


39 


0. 


01956707 


0. 


07540935 


40 


0. 


01957306 


0. 


07517587 
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N 


a 


ai 




0.2 






22 


0. 


01976614 


0. 


05375917 


0. 


08754307 






23 


0. 


01976114 


0. 


05408172 


0. 


08654625 






24 


0. 


01975680 


0. 


05437453 


0. 


08562049 






25 


0. 


01975301 


0. 


05464205 


0. 


08475611 






26 


0. 


01974966 


0. 


05488785 


0. 


08394517 






27 


0. 


01974669 


0. 


05511490 


0. 


08318096 






28 


0. 


01974404 


0. 


05532562 


0. 


08245789 






29 


0. 


01974166 


0. 


05552207 


0. 


08177122 






30 


0. 


01973951 


0. 


05570593 


0. 


08111693 






31 


0. 


01973756 


0. 


05587867 


0. 


08049156 






32 


. 


01973578 


0. 


05604149 


. 


07989215 






33 


0. 


01973415 


0. 


05619545 


0. 


07931614 






34 


0. 


01973265 


0. 


05634145 


0. 


07876131 






35 


0. 


01973127 


0. 


05648027 


0. 


07822571 






36 


0. 


01972999 


0. 


05661259 


0. 


07770765 






37 


0. 


01972881 


0. 


05673898 


0. 


07720563 






38 


0. 


01972770 


0. 


05685998 


0. 


07671833 






39 


0. 


01972667 


0. 


05697603 


0. 


07624459 






40 


0. 


01972571 


0. 


05708754 


0. 


07578337 






TV 


a 


ai 




a 2 




0, 3 




23 


0. 


01970992 


0. 


06084150 


0. 


04428770 


0. 


07428718 


24 


0. 


01971020 


0. 


06080458 


0. 


04451851 


0. 


07388144 


25 


0. 


01971034 


0. 


06078451 


0. 


04464682 


0. 


07365079 


26 


0. 


01971039 


0. 


06077784 


0. 


04469036 


0. 


07357085 


27 


0. 


01971036 


0. 


06078284 


0. 


04465705 


0. 


07363327 


28 


0. 


01971026 


0. 


06079722 


0. 


04455931 


0. 


07382009 


29 


0. 


01971013 


0. 


06081948 


0. 


04440510 


0. 


07412046 


30 


0. 


01970995 


0. 


06084830 


0. 


04420187 


0. 


07452360 


31 


0. 


01970976 


0. 


06088254 


0. 


04395605 


0. 


07501986 


32 


0. 


01970954 


0. 


06092130 


0. 


04367303 


0. 


07560099 


33 


0. 


01970931 


0. 


06096380 


0. 


04335757 


0. 


07625946 


34 


0. 


01970908 


0. 


06100940 


0. 


04301375 


0. 


07698862 


35 


0. 


01970883 


0. 


06105753 


0. 


04264515 


0. 


07778247 


36 


0. 


01970859 


0. 


06110773 


0. 


04225488 


0. 


07863562 


37 


0. 


01970834 


0. 


06115962 


0. 


04184568 


0. 


07954324 


38 


0. 


01970809 


0. 


06121284 


0. 


04141995 


0. 


08050094 


39 


0. 


01970785 


0. 


06126712 


0. 


04097980 


0. 


08150474 


40 


0. 


01970761 


0. 


06132220 


0. 


04052709 


0. 


08255105 
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N 


a 




ai 




0.2 


a 3 




04 




24 


0. 


01971252 


0. 


06028314 


0. 


04946255 


0. 


05630611 


O.i 


32220967 


25 


0. 


01971161 


0. 


06048759 


0. 


04752405 


0. 


06319720 


O.i 


31350151 


26 


0. 


01971081 


0. 


06067477 


0. 


04571025 


0. 


06978710 


O.i 


30499021 


27 


0. 


01971004 


0. 


06086343 


0. 


04384348 


0. 


07671271 


-0 


.00414364 


28 


0. 


01970935 


0. 


06103877 


0. 


04207338 


0. 


08341263 


-0 


.01315896 


29 


0. 


01970871 


0. 


06120830 


0. 


04032868 


0. 


09014499 


-0 


.02239437 


30 


0. 


01970812 


0. 


06137074 


0. 


03862565 


0. 


09683963 


-0 


.03175015 


31 


0. 


01970757 


0. 


06152632 


0. 


03696505 


0. 


10348538 


-0 


.04120539 


32 


0. 


01970707 


0. 


06167592 


0. 


03534052 


0. 


11010004 


-0 


.05078048 


33 


0. 


01970660 


0. 


06181962 


0. 


03375373 


0. 


11666979 


-0 


.06045075 


34 


0. 


01970616 


0. 


06195787 


0. 


03220232 


0. 


12319772 


-0 


.07021604 


35 


0. 


01970575 


0. 


06209097 


0. 


03068515 


0. 


12968223 


-0 


.08006943 


36 


0. 


01970536 


0. 


06221921 


0. 


02920108 


0. 


13612229 


-0 


.09000497 


37 


0. 


01970501 


0. 


06234286 


0. 


02774880 


0. 


14251786 


-0 


.10001832 


38 


0. 


01970467 


0. 


06246219 


0. 


02632719 


0. 


14886860 


-0 


.11010477 


39 


0. 


01970435 


0. 


06257743 


0. 


02493510 


0. 


15517452 


-0 


.12026033 


40 


0. 


01970406 


0. 


06268879 


0. 


02357145 


0. 


16143576 


-0 


.13048133 



We observe in these fits the following behaviour: As N increases, the estimates of 
the amplitudes ai appear to be converging in each set, until we reach the set with five 
asymptotic coefficients (ao, ...,04). In this latter fit, we see that the estimate of 04 
appears to be diverging. Further, the estimates of a±, a<i and 03 have deteriorating 
apparent convergence as we go from a four-term to a five-term fit. By contrast, going 
from a two-term to a three-term fit, and from a three-term to a four-term fit, improved 
the apparent convergence of the amplitude sequences. Thus we reject the five-term fit, 
and base our estimates on the four-term fit. 

Finally, the estimated error from a series of length N is taken to be the appropriately 
scaled difference between the values obtained from this series and those obtained from the 
series of length N — 10. This difference is scaled by a factor dependent on the difference 
between the exponent in question and the first omitted exponent. The scaling factor 
follows from our assumption that the error is given principally by the first neglected 
term, c/N Ak+1 (or similarly with A^J_ 1 ). Hence, if the actual value of the coefficient in 
question is Oj and the two estimates are a\ N ^ and a\ N ~ 10 \ we expect that 

_ 

„ (JV-10) 
(N - 10) A < 
Simple algebra then yields 

(AO 



Therefore, aj is estimated by a K f with error quoted as 

2\af > - af~ »>|/[( 1] . (3.6) 



a,- c 
— — H 

ai 



+ 



(N - 10) A » (N - 10) A *+i 



(3.4a) 
(3.4b) 



(AO (JV-10) 



JV ^A fc+1 -A i _ j 



\ JV-10/ 



(3.5) 
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The factor of 2 is included to make our errors more conservative. More adventurous 
readers may choose to reduce this factor. 

In the above example, the first omitted term is 0(iV~ 5 / 2 ). The difference in the 
estimate of ao from the N = 40 series with four asymptotic coefficients (0.019707611) 
and that from a N = 30 series (0.019709959) is 2.348 xl0~ 6 . Thus the error is quoted 
as 2 x 2.348 x 10~ 6 /((40/30) 2 - 5 - 1) = 4.5 x 10~ 6 . Our amplitude estimate is then 
a = 0.019708 ± 0.000005. Similarly, ai = 0.0613 ± 0.0018, where the error is given by 
2 x 0.0004739/((40/30) L5 - 1). Likewise, a 2 = 0.04 ± 0.02, and a 3 = 0.08 ± 0.1. 

3.3 SAW counts 

In this subsection we discuss the analysis of the newly extended series for SAW counts 
on the square and triangular lattices. Here we give only a brief analysis, as fuller details 
will be published elsewhere [65], along with a discussion of the series derivation. 

Let us begin with the triangular lattice. We first analysed the extended SAP series 
using biased differential approximants using the known exponent a — 1/2 (see [19] for 
the method). We obtained 

Xc = 1/fj, = 0.240 917 574 ±0.000 000 004 , (3.7) 

which we will use in subsequent analyses. 

We also performed a similar analysis using the extended SAW series, biasing the 
estimate with the known exponent 7 = 43/32. We obtained the estimate x c — 1//J, — 
0.240 917 579 ± 0.000 000 008, in agreement with the SAP result but less precise. 

Using the estimate (3.7) of x c , we proceeded as described in Section 3.2 to fit the 
series coefficients to various asymptotic forms. For triangular-lattice SAWs, we expect, 
based on earlier investigations of the corresponding square-lattice series [17], that 

c N ~ ^N 11 ' 32 ^ + di/N + a 2 /N 3 ' 2 + a 3 /N 2 + a A /N^ 2 + a 5 /N 3 + •••]. (3.8) 

However, as discussed in Section 2.2, renormalization-group theory predicts an additional 
"energy-like" term arising from the mixing between nonlinear scaling fields. For the two- 
dimensional SAW (a = 1/2), incorporating this term requires that (3.8) be modified to 
read 

c N ~ pL N N n / 32 [a + ai /N + a 2 /N 3 / 2 + S /iV 59/32 + a 3 /N 2 + a,/N 5 / 2 + • • •] (3.9) 
[cf. (2.26)]. 

In our analysis, we tried both the asymptotic forms (3.8) and (3.9). Since the ex- 
ponents associated with amplitudes 5o and 03 are numerically close (1.84375 and 2, 
respectively), we expect that it will be very difficult to distinguish numerically between 
the Ansatze (3.8) and (3.9). This is indeed the case. Under the assumption (3.8), 
we found that the sequences corresponding to the amplitudes are well converged up to 
k = 4 for triangular-lattice SAW and we estimate a = 1.183966(1), a\ = 0.5960(4), 
a 2 = —0.274(6), a 3 = —0.14(4), and a 4 = 0.09(10). Our errors are calculated as de- 
scribed in Section 3.2 and are given, in parentheses, as the uncertainty in the last quoted 
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digit(s). Under the alternative Ansatz (3.9), we find that the fit is neither better nor 
worse. We observed that the sequences of estimates of the corresponding amplitudes a 
and a 3 appear to be correlated: they are monotonically increasing in magnitude but are 
of opposite sign, the sum a + a 3 being almost constant. 

To investigate this point further, we constructed a test series, with known asymptotic 
behaviour, similar to that in (3.9), namely 

d N = N u/32 [1 + 1 /N + . 7 /N 3/2 + 1.25 /iV 59/32 + 3/iV 2 - 4 /N 5/2 + 5/iV 3 - 6 /N 7/2 } + (0.5) N . 

(3.10) 

The last term is included to incorporate the fact that there are other singularities in the 
complex plane, beyond x c , which will make an exponentially decaying contribution to 
the asymptotics. We generated the first 1000 terms of this sequence and analysed them 
as above, including either a term jV _59//32 or a term N~ 2 or both. The analyses using 
either one of these two terms behaved similarly. The analysis using both terms gave 
inferior estimates of the first three amplitudes, and the wrong sign for the amplitude 
of the N~ 2 term, when N < 240. Only beyond this point does the analysis using both 
terms give superior estimates of the first three amplitudes, along with the right sign for 
the N~ 2 term (the two issues clearly go together). We conclude that using series of the 
length available to us (N < 40), it is unfeasible to determine whether a term TV -59 / 32 is 
present or absent. 

In conclusion, our analysis is unable to resolve the question of whether the "energy- 
like" term jV~ 59//32 is present or not. Therefore, for the subsequent analysis of the metric 
quantities, reported in the next subsection, we have assumed for simplicity the absence 
of this term, and just assumed the asymptotic form (2.24) with one correction-to-scaling 
exponent Ai = 3/2. 

On the square lattice, the situation is complicated by the presence of an "antiferro- 
magnetic" singularity at x = — l//i. From (2.27) ff. we recall that the asymptotic form 
of the coefficients given in (3.8) and (3.9) is modified by the additional term 

(-1) N fi N N- 3 / 2 [d + h/N + d 2 /N 2 + d 3 /N 3 + •••]. (3.11) 

Analysing the square-lattice SAW data using as our estimate of x c the positive real root 
of 

581x 4 + 7x 2 - 13 = 0, (3.12) 

which is a useful mnemonic for the current best estimate x 2 = 0.143 680 629 27(1) [19], 
we obtain a very convincing fit with k = 3 and m — 2, enabling the following amplitude 
estimates to be made: a = 1.1770425(7), a ± = 0.5501(2), a 2 = -0.1402(3) and a 3 = 
-0.12(2) for the "ferromagnetic" amplitudes, and d = -0.189848(3), d x = 0.17473(9) 
and g?2 = —1.51(1) for the "antiferromagnetic" amplitudes. This is in close agreement 
with the earlier estimates (1.2) based on slightly shorter series; here we have obtained 
a slight improvement in the precision of the estimates for the leading amplitudes a , d 
and d\. 

If instead we assume the asymptotic form (3.9) for the "ferromagnetic" term, we find 
that estimates of a 3 are small (less than 0.03 in magnitude) and tending toward zero. 
Estimates of a 4 are tending toward the estimate for this term obtained in the absence of 
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the additional term with amplitude 03. Once again, we are unable to distinguish (3.8) 
from (3.9). 

3.4 Metric quantities (R 2 e ), (R 2 g ) and (R 2 m ) 

In this subsection we shall analyse the metric quantities (R 2 )n, (Rq)n and (i?^)jv, any 
one of which we shall generically denote by (R 2 )n- As discussed in Section 2.2 [cf. 
(2.24)/(2.27)], their asymptotic behaviour is expected in the first instance to be 

(R 2 ) N ~ N 2u [a Q + a±/N + a 2 /iV Al + a 3 /N 2 + ■ ■ ■] (3.13) 

for the triangular lattice, and 

(R 2 ) N ~ iV 2 "[a + ai/iV + a 2 /iV Al + a 3 /iV 2 + • • •] 

+ (-1)^^(60 + &i/iV A ^ + & 2 /iV + • ■ ■] (3.14) 

for the square lattice. As mentioned earlier, there is overwhelming numerical evi- 
dence [17,37] that the leading exponent 2v equals 3/2 exactly, as predicted by Coulomb- 
gas arguments [4,5]; we shall henceforth take this fact for granted. We also expect 
Ai = 3/2, as predicted by Nienhuis [4] and confirmed numerically for the SAW counts. 
Furthermore, for the square lattice we predict q = —11/32 [cf. (2.33)]. Finally, in 
Section 2.2 we pointed out that renormalization-group theory predicts an additional 
"energy-like" term in the susceptibility [cf. (3.9)], though alas we were unable to dis- 
tinguish it numerically (see Section 3.3); it is reasonable to guess that there may be a 
corresponding term also in the series for the unnormalized second moments C]y{R 2 )iy. 
Whether or not the latter term is present, the existence of an "energy-like" term in c N 
will induce in (R 2 )n additional correction-to-scaling terms jV~ 59//32 , jV~ 91//32 , . . . beyond 
those included in (3.13)/(3.14). 

In addition to the normalized metric quantities (R 2 )n, we also studied the corre- 
sponding unnormalized quantities cjy(R 2 )]y, whose expected asymptotic form is 

c N (R 2 ) N ~ ^W+^t&o + bjN + 6 2 /iV Al + b 3 /N 2 + •••], (3.15) 

with appropriate additional antiferromagnetic terms (2.31) when analysing the square- 
lattice data. The latter quantities have the disadvantage that the analysis depends 
sensitively on an input estimate of //; but, for loose-packed lattices and for (R 2 )n only, 
they have the advantage that the effect of the antiferromagnetic singularity is weaker. 
To see this, compare (2.30)-(2.33): the antiferromagnetic contribution in cn(R 2 )n is 
relatively weaker than that in cn] but the antiferromagnetic contribution in (R 2 )n is 
dominated by that in cat. Therefore, the antiferromagnetic contribution is relatively 
weaker in cn(R 2 )n than in (R 2 )n- 

Our method of analysis is based on directly fitting (R 2 )n and cn(R 2 )n to the as- 
sumed asymptotic form (3.13)/(3.14)/(3.15), as described in Section 3.2. The values 
of the exponents u, q, Ai and Af F are assumed, and the appropriate system of linear 
equations is solved to give estimates of the amplitudes {ai} and In applying the 

method to metric quantities (see, for example, the table in Section 3.2), the fit to the 
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leading amplitude is rather stable, that to the first analytic correction term is moderately 
stable, while the fit to the amplitude of the assumed correction-to-scaling term N~ 3 ^ 2 
converges less impressively for both the normalized and unnormalized metric quantities. 
As already noted, adding further terms in the assumed asymptotic form beyond the first 
initially improved convergence, but this improvement is not sustained. That is to say, 
allowing more than between two and five terms in the assumed asymptotic form led to 
an apparently deteriorating fit. 

3.4.1 Triangular lattice 

With our 40-term triangular-lattice series we found that we could fit to a , a±, a 2 and 
sometimes a 3 . For the normalized and unnormalized metric quantities, the estimate of 
a could usually be made to four-digit precision, while the estimate of a\ could be made 
only to one or two significant digits, and the estimate of a>2 is accurate only to at best 
one significant digit. For a 3 the error is comparable to or greater than the estimate. 

We have applied this analysis method to the triangular-lattice data, using a 40-term 
series for all metric quantities. Because the triangular lattice is close-packed, there is 
only one singularity on the circle of convergence, which makes the analysis simpler than 
for the square lattice [compare (3.13) to (3.14)]. 

The tables of estimates for the metric quantities obtained according to the procedure 
described in Section 3.2 are shown in Table 7. The results are: 

(R 2 e ) N ~ N 3/2 [0.71 174(32) + 0.95(12)/JV - 2.6(1.6)/iV 3/2 + 3(7)/iV 2 + 0(l/iV 5/2 )] 

(3.16) 

(R 2 g ) N ~ N 3/2 [0.09989(4) + 0.056(16)/iV + 0.3(2)/iV 3/2 - 0.2(1. 0)/iV 2 + 0(l/iV 5/2 )] 

(3.17) 

(R 2 J N ~ iV 3/2 [0.3133(4) + 0.24(12)/A^-0.2(1.0)/A^ 3/2 + O(l/iV 2 )]. (3.18) 

These were obtained with k — 3, k — 3, and k = 2 respectively. Unfortunately the 
uncertainties in the coefficients of the 0(l/iV 3//2 ) are so great as to be comparable to (or, 
in the case of (-R 2 ^), larger than) the coefficient itself. Further, the analysis of the {R 2 )n 
series violates the convergence criterion we have set, in that the coefficient of a\ differs 
by nearly a factor of 3 in going from a two-term fit (k = 1) to a three-term fit (k = 2). 
We have nevertheless presented results for k — 3. Our justification for this is twofold. 
Firstly, the estimate of a x stabilises if we then go to a four-term fit. Secondly, as we have 
already seen, the data for the essentially equivalent series (N+l) 2 c N (R 2 ) N /Q supports a 
four-term fit. For the reader unconvinced by these arguments, the corresponding analysis 
with k = (a one-term fit, as would be justified by strict adherence to the convergence 
criteria we have set), gives ao = 0.106 ± 0.012. 

We can do somewhat better from a similar analysis of the unnormalized metric 
quantities, using the estimate \i = 4.15079723 from (3.7), which gave 

c N (R 2 e )N/Q ~ /i JV A^ 59/32 [0.14045(6) + 0.256(26)/A^- 0.53(32)/A^ 3/2 

+0.6(1.6)/iV 2 + O(l/iV 5/2 )] (3.19) 
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(N + l) 2 c N (R 2 g } N /Q ~ fi N N 123/32 [0.019708(5) + 0.0613(18)/A^ + 0.04(2)/iV 3/2 

+0.08(10)/iV 2 + 0(l/iV 5/2 )] (3.20) 
(N + l)c N (R 2 J N /Q ~ //^ 91/32 [0.06183(12) + 0.136(28) /iV - 0.02(24) /N 3/2 + 0(1/N 2 )} 

(3.21) 

These estimates were also obtained with k = 3, k = 3 and k = 2 respectively. From 
these we can estimate the amplitudes of the metric quantities by dividing through by 
the asymptotic form (3.8) for c/v, and accounting for the factor of 6 and appropriate 
factors of (N +1). In this way we obtain a second set of amplitude estimates, 

(R 2 e } N ~ iV 3/2 [0.71 176(15) + 0.94(7) /N - 2.5(8)/iV 3/2 + 3(3)/iV 2 + 0(l/iV 5/2 )] 

(3.22) 

{R 2 g ) N ~ N 3/2 [0.09987(2) + 0.061(5)/iV + 0.23(5)/N 3/2 - 0.5(5)/N 2 + 0{l/N b/2 )] 

(3.23) 

(R 2 m ) N ~ iV 3/2 [0.3133(3) + 0.22(7)/iV-0.03(60)/A^ 3/2 + O(l/iV 2 )]. (3.24) 

These differ from the directly measured amplitudes only within the quoted errors for 
each amplitude, consistent with our claimed errors. 

One immediate observation is that for (R 2 )n and (R^n the correction-to-scaling 
amplitudes corresponding to the 1/N and 1/iV 3 / 2 terms are of opposite sign, while for 
(R 2 ) A? they are of the the same sign (though the errors associated with the estimates of 
the amplitude of the 1/N 3 / 2 term are rather large). Note too that for both (R 2 )n and 
(R 2 )n the amplitude of the 1/iV 3 / 2 term is larger (in magnitude) than the amplitudes 
of both the leading term and first analytic correction; for (R^), by contrast, the error 
in the 1/N 3 / 2 term is too great to comment on the relative size of this term. 

As a consequence, the "effective" exponent A cff based on fitting to a given range of 
A" behaves differently as a function of A" for the different observables. For (i? 2 )jv, A e fr 
lies between 1 and 3/2 for all A", and decreases monotonically to 1 as A" — > oo. For 
{R 2 )n and (R^n, by contrast, A cff starts above 3/2 for very small N, then increases 
monotonically, reaching +oo at some finite value of A^ (here ~ 15); then it jumps to 
— oo, after which it continues to increase monotonically, tending asymptotically to 1 as 
A^ — > oo. These observations are in accordance with previous studies. Most studies 
of (R 2 )n resulted in estimates for A of 0.65 [8,10,13,14,23], while most studies of 
(R 2 )n resulted in estimates for A in the range 1.05-1.2 [13,14,20,23]. (There have 
been few previous studies of (R 2 m )N.) This is clearly a source — indeed, probably the 
major source — of the long-standing difficulty in the analysis of these quantities for the 
correction-to-scaling term. 

The amplitude ratios A and B [defined as the A^ — > oo limit of (2.8)] follow imme- 
diately from (3.22)-(3.24) as A = 0.14031 and B = 0.4398, and from (3.16)-(3.18) as 
A = 0.14033 and B = 0.4402. 

We also estimated the ratios A and B by direct extrapolation of the appropriate 
coefficient quotients, using the following method [41]: Given a sequence {a n } defined for 
n > 1, which is known or assumed to converge to a limit with corrections of the form 
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a n ~ a 00 {\+b/n + . . .), we first construct a new sequence {h n } defined by h n = Yim=i a m- 
Then the generating function J2h n x n ~ (1 — a^x)"^ 1 ^. We obtain estimates of the 
required limit and parameter b by analysing this generating function by the standard 
method of differential approximants. (The value of the parameter b can also be obtained 
numerically from the amplitude estimates given in (3.22)-(3.24) above.) In this way, we 
obtain the estimates 

A = 0.140296(6) (3.25) 
B = 0.439649(9) . (3.26) 

3.4.2 Square lattice 

Let us now consider the square-lattice data. We first analysed the three metric quan- 
tities (-Kg) at, (R 2 )n and (-R^)at by a method similar to that leading to (3.16)-(3.18), 
but including the contribution of the antiferromagnetic singularity. We imposed the 
exponent values v = 3/4, Ai = 3/2, q = —11/32 and Af F = 1; the justification for 
these choices has already been given above. The sequences of amplitude estimates are 
shown in Tables 8, 9 and 10. In this way, we obtain the following results: As discussed 
in Section (3.2), some experimentation was needed to determine the maximum values of 
the parameters m and k in (3.2). The results given below are given by m = 4, k = 2 for 
(RI)n, by m = 2, k = 1 for (FQn, and by m = 3, k = 1 for {R 2 m )N- We find 

(-Re) n ~ N 3/2 [0.77124(5) + 1.159(38)/7V - 3.13(74)/iV 3/2 + 6(6)/iV 2 - 6(24) /N 5/2 

+ 0.4(4.0)/iV 3 + O(l/iV 7/2 )] + (-l) Jv iV- 11/32 [0.12451(17) - 0.027(24)/iV + 0(l/iV 2 )] 

(3.27) 

(R 2 g ) N ~ iV 3/2 [0. 108230(1) + 0.1019(l)/iV + 0.1082(4) /iV 3/2 - 0(l/iV 2 )] 

+ (-l) N N~ n/32 [0.008364(19) + 0.0031(21)/JV + 0(l/iV 2 )] (3.28) 

(/Oat ~ A^ 3/2 [0.33913(8) + 0.426(17) /N - l.l(l.l)/iV 3/2 + 2(4)/A^ 2 + 0(l/iV 5/2 )] 

+ (-l) Ar iV- 11/32 [0.03652(ll) + 0.015(12)/7V + 0(l/iV 2 )] . (3.29) 

A similar analysis of the unnormalized quantities, using the estimate 
pi = 2.63815853034174086843 from (3.12), was made. The results below are given by 
m = 4, k = 2 for the first two quantities, and m = 3, k = 1 for the third. For the first 
two quantities we have not given our estimate of 04 as the associated error is significantly 
bigger than the estimate. We find 

c N (R 2 e ) N /A ~ iV 59/32 [0.226945(14) + 0.4471(ll)/iV - 0.95(22)/iV 3/2 + 2(2)/iV 2 + 0(l/iV 5/2 )] 

+ (-1) N N~ 3/2 [0.019098(1) + 0.0415(41)/JV - 0.08(45)/iV 2 + 0(l/iV 3 )] (3.30) 

c n {R 2 9 )n ~ A^ 59/32 [0.127388(31) + 0.181(17) /N + 0.10(26)/A^ 3/2 + 0.1(1. 0)/A^ 2 + 0(l/iV 5/2 )] 
+ (-1)^-0.010688(15) + 0.0047(17) /N - 0.20(5)/iV 2 + 0(1/A^ 3 )] (3.31) 

cn(R 2 Jn ~ A^ 59/32 [0.39917(ll) + 0.686(44)/iV - 1.3(6)/iV 3/2 + 2(2)/iV 2 + 0(l/iV 5/2 )] 

+ (-1) JV [-0.021383(45) + 0.028(5)/A^ + 0(l/iV 2 )] . (3.32) 
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From these we can estimate the amplitudes of the normalized metric quantities by di- 
viding through by the asymptotic form (2.30) for cjy. In this way we obtain a second 
set of amplitude estimates, 

(R 2 e ) N ~ N 3/2 [0.77124(5) + 1.1593(38)/iV - 3.13(75)/iV 3/2 + 6(8)/iV 2 + 0(l/iV 5/2 )] 
+ (-l) iY ^ 11/32 [0.12439(4) - 0.0144(9)/A^ + O(l/iV 2 )] (3.33) 

(R 2 g ) N ~ iV 3/2 [0.108227(50) + 0.103(17)/iV+ 0.098(160)/A^ 3/2 -O(l/iV 2 )] 

+ (-l)^^- 11 / 32 [0.008376(20) + 0.0006(30)/iV + 0(l/iV 2 )] (3.34) 

(R 2 m ) N ~ iV 3/2 [0.33913(8) + 0.424(19)/iV-l(l)/iV 3/2 + O(l/iV 2 )] 

+ (-l) N N- 11/32 [0.03654(7) + 0.027(8)/iV + 0(l/iV 2 )] . (3.35) 

These differ from the directly analysed amplitudes only in the last quoted digits for all 
but the least significant amplitudes, and are consistent with our quoted errors in all 
cases. In the notation of (2.34), from (3.33) we have a'^/dl = 6.200(3), while from the 
amplitudes quoted below (3.12) we have do/do = 6.1999(1), in complete agreement. 

The amplitude ratios A and B follow immediately from (3.27)-(3.29) as A = 0.14033 
and B = 0.43971. From (3.33)-(3.35) we obtain the almost identical values, A = 0.14033 
and B = 0.43972. 

We also analysed these amplitude ratios directly, using the same method as discussed 
above for the analysis of the triangular-lattice data. We obtained 

A = 0.140299(6) (3.36) 
B = 0.439647(6) . (3.37) 

Comparison with the corresponding estimates (3.25)-(3.26) for the triangular lattice is 
entirely consistent with the belief that these ratios are lattice-independent [25]. 

These amplitude ratios are also consistent with the CSCPS relation (2.10): using our 
best estimates (3.36)-(3.37), we find F = lim^oc^ = -0.000024 ± 0.000025 for the 
square lattice; and using (3.25)-(3.26), we find F = lim^oc F N = -0.000036 ±0.000036 
for the triangular lattice. A direct analysis of the sequence {Fn} was also undertaken, 
but that sequence was found difficult to extrapolate; and our estimate of the limit, while 
entirely consistent with zero, was a factor of 10 less precise than the one just quoted. 

If it is in fact true (as certainly seems to be the case) that Fn — > as N — > oo, 
then it is of interest to investigate the rate at which Fn — > 0. If Fn oc N~ s , then 
Jn = Fn{R 2 )n oc jV 3 / 2-5 . For both square and triangular lattices, we find that 5 = 3/2, 
i.e. that /n approaches a nonzero constant / = hin^^oo f N &sN^ oo. This behaviour 
was initially surprising to us, because it implies that there is no analytic correction-to- 
scaling term 1/N in Jn, even though such terms are manifestly present in each of the 
three individual metric quantities (-R 2 ). Moreover, this remarkable result appears to hold 
for both lattices. 6 However, we were subsequently able to provide a renormalization- 

6 Indeed, as shown in [41] it appears to hold even in the case of interacting SAWs within the good- 
solvent regime (i.e., above the theta temperature). Of course, the limiting constant / depends on the 
interaction. For repulsive nearest-neighbour interactions, / increases from 0.78 to an asymptotic value 
of about 1.6 as the repulsion gets very strong. 
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group argument for this cancellation (see Section 2.2 above). Our estimates of the 
amplitude are / = 0.78 ± 0.03 [64] for the square lattice and / = 0.96 ± 0.04 for 
the triangular lattice. These estimates are based on extrapolation of the sequences f N 
using a variety of extrapolation algorithms, including Levin's u transform, Brezinski's 6 
algorithm, Neville-Aitken extrapolation and Wynn's e algorithm. Details of these and 
other algorithms, as well as programs for their implementation, can be found in [70]. 

As noted in Section 2.2, the observation that 5 = 3/2 implies the constraint (2.44) on 
the subdominant amplitudes arising in (2.40)-(2.42). Our series estimates (3.22)-(3.24) 
and (3.27)-(3.29) are consistent with this prediction, as are our Monte Carlo estimates 
(4.3)-(4.11). Furthermore, our series and Monte Carlo estimates of / are consistent with 
the relation (2.45); but the associated error bars are very large, so this is not a stringent 
test. 

Finally, we note the fact that 5 = 3/2 is another indicator that the correction-to- 
scaling exponent is indeed 3/2. If it were less than this, then the leading non- analytic 
correction-to-scaling term would have to cancel miraculously (as the 1/N term does) in 
the combination (2.11) for f^. This seems a priori unlikely. 

3.5 Euclidean-invariant moments of the distribution function 

We have also analysed the series for rotationally-invariant and non-rotationally-invariant 
moments of the endpoint distribution function, given in Tables 3-6, using methods 
similar to those just described for the analysis of (R 2 )- Let us start with the rotationally- 
invariant moments {r 2k )^ 1 for which we expect an asymptotic behaviour of the form 

{r 2k ) N ~ N 2kv [c 0>k + c 1>k /N + c 2>k /N Al +c 3 , k /N 2 + ...] (3.38) 

for the triangular lattice, and 

(r 2k ) N ~ N 2k »[c , k + c hk /N + c 2 , fc /iV Al + c 3 , k /N 2 + . . .] 

+ (-1) N N«* [d 0) * + d hk /N A > F + d 2 , k /N +...} (3.39) 

for the square lattice. In Section 2.2 we gave arguments predicting that q k = 2kv + a — 
1 — 7 = 3k/2 — 59/32. Furthermore, it is reasonable to expect A^ F = 1 as was already 
observed for the SAW counts [17] and for the metric quantities (-R 2 ). Our numerical 
results are consistent with these predictions. 

We began by analysing the moments (r 2k )^ using the method of differential approx- 
imants [70,71], with the aim of confirming the predicted leading exponents 2kv = 3k/2 
and q k = 3k/ 2 — 59/32. It is a previously observed feature of the method of differential 
approximants (DA) that its application to the analysis of SAW moment series is less 
accurate than might be expected [71]. For example, DA analysis of a 27-term square- 
lattice (R 2 )n series, biased at a critical point of 1, produces estimates of 2v in the range 
1.495-1.497 [71], whereas an unbiased analysis of a 27-term SAW series on the same lat- 
tice yields exponent estimates of 7 = 1.34364 ± 0.00088, which is rather more accurate, 
as well as more precise. This behaviour is most likely connected to the fact that the 
method of differential approximants tacitly assumes that the function is well approxi- 
mated by a differentially finite (D-finite) function, i.e. the solution of a linear ordinary 
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differential equation with polynomial coefficients [72]. While there is strong evidence 
that neither SAWs nor S APs are D- finite [72] , it nevertheless appears that the SAW and 
SAP counts are well approximated by a D-finite function, while the generating functions 
for SAW and SAP metric properties (such as {R 2 ) n) appear not to be. Evidence for 
this remark includes the telling fact that most of the differential approximants for (R 2 ) 
(for both SAWs and SAPs) are defective [71], which is usually a signal that the function 
being approximated is not of the type tacitly assumed by the analysis. For this reason, 
our DA analysis gives only moderately accurate estimates of the leading exponents, both 
ferromagnetic and antiferromagnetic, but no reliable information as to the value of the 
subleading exponents. 

Our DA analysis confirmed the expected leading behaviour (r 2k )]y oc jV 3fc / 2 , the 
exponents being identified as 1.4997(5), 2.998(6), 4.496(9), 5.996(12), and 7.496(12) for 
k = 1,2,3,4,5 on the square lattice, and 1.4997(5), 2.996(7), 4.495(9), 5.997(12), and 
7.500(10) for k — 1,2, 3, 4, 5 on the triangular lattice. DA analysis also gave reasonable 
estimates of the leading antiferromagnetic exponent on the square lattice: we found 
qk = 3k/2 — 59/32 to an accuracy of approximately 0.01, 0.04, 0.05, 0.05, 0.06 for 
k = 1,2,3,4,5, respectively. It is possible that higher moments (k > 1) may have non- 
analytic correction-to-scaling terms with exponent Ai < 1.5 which would then be more 
prominent than the leading correction-to-scaling term of the second-moment quantities 
(R 2 )- A more prominent such singularity would also explain the relatively poor accuracy 
of the DA analysis. We allow for this possibility in our analysis, described immediately 
below, but find no evidence for such a term. 

We next proceeded to fit (r 2k ) N to the asymptotic forms (3.38) and (3.39), setting 
v — 3/2 and = 3k/ 2 — 59/32 and investigating the quality of the fit for a variety of 
possible values of Ai and A^ F . Among other things, we considered the possibility that 
Ai < 1, even though no such term is observed in the metric quantities (R 2 )n- 

We first fitted the available series to the above forms with Ai = 1/2. Estimates 
of the associated amplitude were, in all cases, monotonically decreasing toward zero. 
Furthermore, as we increased the number of subdominant terms included in the fit, this 
amplitude approached zero more and more closely. The data are insufficient to judge 
whether the rate of approach to zero increased, but the entries were numerically smaller. 
Also, estimates of the leading amplitude co,k did not display the sort of convergence we 
found in the analysis of {R 2 )n', rather, the convergence deteriorated as we increased the 
number of subdominant terms included in the fit. Both of these observations suggest 
that there is no term C2,k/N Al with Ai pa 1/2. Similar behaviour was observed with 
Ai = 11/16, though we cannot say whether the effect was stronger or weaker. That is 
to say, the analysis is insensitive to this level of exponent change for these series. This 
is consistent with the situation found in the analysis of (R 2 ). We conclude that there is 
no evidence of a correction term C2,k/N Al with A x < 1. 

We then reanalysed the data assuming that the only correction-to-scaling term, other 
than integer powers of 1/N, was that with exponent Ai = 3/2, exactly as found for 
the second-moment series. As for the second-moment series, we retained only analytic 
correction terms to the antiferromagnetic singularity. As we increased the order of the 
fit, the leading amplitudes c 0j fc and d 0jk displayed improved convergence. This is usually 
an indicator that the guessed asymptotic form is correct. The higher-order amplitudes 
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displayed less convincing convergence, but we ascribe this to a lack of adequate data. For 
the second moment, we have a 59-term series, which converged rather well, as can be seen 
from Table 8. But the convergence is much less impressive after only 32 terms — which 
is all we have available for the higher moments. 

Taken together, our results favour the most obvious conjecture, which is that the 
subdominant behaviour is characterised by the same exponent set as is observed for 
(R 2 )n- 

In order to test the conclusion that the leading correction term in all these series 
is the 1/N term, we used the method of differential approximants on a modified series 
obtained by subtracting the estimated leading-order term from the original series: that 
is, we analysed (r 2k )N — co t kN 2ku . We found that the series coefficients behave like 
jy2fci/-i^ cons i s tent with the conclusion that the leading correction term is 1/N and that 
the non-analytic correction-to-scaling term(s), have exponent Ai > 1, consistent with 
our view that A ± — 3/2. 

With the foregoing observations in mind, we obtained the following estimates for 
the corresponding amplitudes for the square-lattice moments (r 2k )N, where we have 
assumed a single correction-to-scaling exponent Ai = 3/2 associated with the ferromag- 
netic singularity, and otherwise only analytic corrections to both the ferromagnetic and 
antiferromagnetic singularities: 



k 


= 2: 


c 0,2 


= 0.860(2), ci, 2 = 


1.9(2), 








do,2 


= 0.139(5), d 1>2 = 


-0.03(2), 


(3.40) 


k 


= 3: 


Co,3 


= 1.184(5), ci, 3 = 


3(1), 








^0,3 


= 0.193(5), d 1>3 = 


-2(1), 


(3.41) 


k 


= 4: 


c 0,4 


= 1.907(10), c M = 


= 2.5(5), 








do,A 


= 0.310(3), d 1>4 = 


-0.46(5), 


(3.42) 


k 


= 5: 


Co,5 


= 3.434(10), ci )5 = 


= -3(1), 








do,5 


= 0.551(3), di >5 = 


-1.6(2). 


(3.43) 



As a check we verify relation (2.34). Our results for cn predict do,k/co,k = 0.161292(2), 
a relation that is well satisfied by our results for co,fc and do,fc- 

We can now provide a direct estimate of the invariant ratios M 2 k,N = (r 2k ) N / (r 2 )% in 
the limit N — > oo. From the above amplitude estimates, we have, for the square lattice, 



M 4i00 


= 1.446(3) 


(3.44) 


M 6iOC 


= 2.581(11) 


(3.45) 


M 8j00 


= 5.391(28) 


(3.46) 


M 10>oo 


= 12.59(4). 


(3.47) 



These estimates agree well with the estimates (2.18)-(2.21) obtained previously [42], but 
are a factor 5-10 less precise. 
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For the triangular lattice, there is of course no "antiferromagnetic" singularity, so that 
the terms corresponding to the amplitudes dj k are absent. We find from the triangular 
lattice data: 



k = 2: c , 2 = 0.7330(9), c lj2 = 1.2(2), c 2 , 2 = -5(1), (3.48) 

k = 3: c ,3 = 0.934(2), c 1>3 = 1.5(5), (3.49) 

k = 4: c ,4 = 1.383(3), c M = 2(1), (3.50) 

k = 5: c ,5 = 2.31(3), c 1>5 = -8.4(5). (3.51) 
From the above amplitude estimates, we have, for the triangular lattice, 

M 4j00 = 1.446(2) (3.52) 

M 6)OC = 2.588(5) (3.53) 

M 8)00 = 5.381(12) (3.54) 

M 10 ,oo = 12.64(15). (3.55) 



These estimates agree well with those found for the square lattice, confirming the ex- 
pected universality. Therefore they are also in agreement with the field-theory estimates 
(2.18)-(2.21), though less precise. 

3.6 Non-Euclidean-invariant moments of the distribution func- 
tion 

In this section we discuss the behaviour of the non-rotationally-invariant moments: 
(r 2k cos A9) N with k = 2,3,4 and (r 8 cos88) N for the square lattice, and (r 2k cosQ6) N 
with k — 3, 4 for the triangular lattice. 

Let us first consider the triangular lattice. We began by analysing the series using the 
method of differential approximants, with the aim of determining the leading exponent. 
For the triangular lattice we write 

(r 2k cos Q9) N ~ jV 2fc "- A " [a , k + a lik /N A + a 2 , k /N + •••]. (3.56) 

Using first- and second-order differential approximants, we found A nr = 3.00 ± 0.10 for 
k = 3 and A nr = 2.95 ± 0.10 for k = 4, from which we conjecture that A nr = Av = 3 
exactly, as predicted in Section 2.1. 

Fitting the triangular-lattice data to the asymptotic form (3.56), we found good 
convergence only if we set the correction-to-scaling exponent to a value A ps 1/2— 
in sharp contrast to situation for the corresponding rotationally invariant moments, 
where we found A = 3/2. To test the conjecture that the leading correction is N^ 1 / 2 , 
we subtracted the estimated leading term ao 

N 2kv-A nl from 

(r 2h cos 66} jv and analysed 

the resulting series. It was found to behave as ai ^,jV" 2fcly ~ Anr ~ a50±a10 , implying that 
A = 0.5 ± 0.10. At this stage, we have no theoretical explanation for this numerical 
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observation. Setting A nr = 3 and A = 1/2 and assuming subsequent half- integer terms in 
the asymptotic expansion (3.56), we obtained the following estimates for the triangular- 
lattice amplitudes: 

k = 3: a , 3 = 1.120(3), a 1>3 = -1.95(5), a 2 , 3 = 1.7(4) (3.57) 

k = 4: a , 4 = 4.05(5), a M = -9(1), a 2 , 4 = 20(4). (3.58) 

For the square lattice, equation (3.56) needs to be modified by the addition of a term 
representing the antiferromagnetic singularity, so we write 

(r 2k cosA9) N ~ N 2k »- A ™[a , k + a hk /N A + a 2 , k /N + ---] 

+ (-l)^[6o lfc + ht/N 1 + b 2 /N 2 + •■■]. (3.59) 

From first- and second-order differential approximants applied to the square-lattice data, 
we found A nr = 1.46 ± 0.03, A nr = 1.45 ± 0.06 and A nr = 1.44 ± 0.09 for k = 2,3,4, 
respectively. From these results we conjecture that A nr = 2v = 3/2 exactly, as predicted 
in Section 2.1. Differential-approximant analysis also gave reasonable estimates of the 
the leading antiferromagnetic exponent: we found ip = 2kv — A nr — 3+7 = 2kv — 101/32, 
accurate to ±0.013 for k = 2, ±0.05 for k — 3, and ±0.15 for k — 4. This expression 
for ip is different from the one which is the natural generalization of the result for the 
Euclidean-invariant moments, 2kv — A nr ± a — 1 — 7 = 2kv — 107/32, which is excluded 
from the analysis: the difference is 6/32 = 0.1875, much larger than the errors. 

Fitting the data to the above asymptotic form (3.59) with A nr = 3/2 and ip = 
2kv — 101/32, and assuming only analytic corrections to scaling at the antiferromag- 
netic critical point, as found for all the other series, we again found good convergence 
only if we set the ferromagnetic correction-to-scaling exponent A to approximately 1/2, 
just as was found in the analysis of the triangular-lattice data. As in the triangular- 
lattice analysis, we verified this conjecture by subtracting the estimated leading term 
a o,fc-^ 2feiy ~ Anr from (r 2k cos 49) jv and analysing the resulting series, which was found to 
behave as ai ^jV 2fciy ~ Anr ~ a497±a005 . This is strong support for an iV -1 / 2 correction. 

Setting A nr = 3/2, ifj = 2kv — 101/32 and A = 1/2 and assuming subsequent half- 
integer terms in the asymptotic expansion of the ferromagnetic singularity (3.59), we 
obtained the following estimates for the square-lattice amplitudes: 



k = 2: 


a 0,2 


= 1.148(6), ai, 2 = -1.70(5), a 2 , 2 = 2.7(3), 






bo,2 


= 0.060(5), 6i, 2 = 0.6(2), 


(3.60) 


k = 3: 




= 3.200(10), ai, 3 = -6.25(10), a 2 , 3 = 12.0(5), 






bo,3 


= 0.175(10), 61,3 = 1.2(4), 


(3.61) 


k = 4: 


O0,4 


= 8.90(10), a M = -20(2), a 2j4 = 50(8), 






b ,4 


= 0.47(5), b 1A = 3(1). 


(3.62) 



Finally, we analyzed the series (r 8 cos 89) jv on the square lattice. A differential- 
approximant analysis gave (r 8 cos89) N oc jV 3 - 07±ai . We conjecture that the exponent is 
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exactly 3, consistent with the behaviour (r 8 cos 80) n oc jV 8iy ~ Anr ' 8 with A nr;8 = 3. This is 
exact at the mean-field level and also for the two-dimensional Ising model. In the same 
spirit that we previously conjectured that A nr = 2u, we now conjecture that A nr 8 = 4za 
The antiferromagnetic exponent in (r 8 cos86)^ was estimated to be ip = 1.35 ± 0.05, 
which, by analogy with the antiferromagnetic exponents for (r 2fc cos40)jv, we conjecture 
is exactly Su — A nrj8 — 3 + 7 = 43/32. We found the subsequent analysis consistent with 
only analytic corrections at the antiferromagnetic critical point. At the ferromagnetic 
critical point, the data were again consistent with a leading N~ l l 2 correction-to-scaling 
term. In an identical notation to that used above, we find the amplitudes to be: 

a = 135(2), ai = -540(10), a 2 = 1440(50), 

b = 6.5(1), 6i = -11(3). (3.63) 

4 Monte Carlo analysis 

4.1 Summary of our data 

We have also generated extensive Monte Carlo data, using the pivot algorithm [37,73-75], 
for SAWs on the square lattice at selected values of N in the range 40 < N < 4000, 
measuring the observables {R 2 ), (R 2 g ) and (R^)- The results are collected in Tables 11 
and 12. Unfortunately, in some runs we did not measure all observables: in particular, for 
larger values of N, the statistics available for (R^) are much smaller than for the other 
radii. The statistics range from 10 9 to 10 10 pivot iterations per point, or approximately 
10 7 to 10 9 times the integrated autocorrelation time r in t,A for these observables. 7 

We begin by analyzing the Monte Carlo data in an unbiased way, in order to ex- 
tract the leading amplitudes and the correction-to-scaling exponents and amplitudes 
(Section 4.2). Then we compare the Monte Carlo data, which lie at relatively high 
iV but are afflicted by statistical errors, with the formulae (3.27)-(3.29) obtained by 
extrapolating the series data, which lie at much smaller iV but are exact (Section 4.3). 

4.2 Data analysis 

In order to determine the leading amplitudes and the correction-to-scaling exponent (s) 
and amplitude(s), we have analysed the three quantities (R 2 ), (R 2 ) and (R^)- We first 
tried nonlinear least-squares fits of the form 8 

{R 2 ) N /N 2 » = a + b(N/N o r A , (4.1) 

7 High-precision Monte Carlo data for {R^)n have been kindly provided by Peter Grassberger. His 
results have been merged with ours and appear in Tables 11 and 12. 

8 Note that we have rescaled the length JVbya fixed parameter N that has always been taken equal 
to N = 750. The purpose of this rescaling is to minimize the covariance between the estimates of b and 
A. (The optimal choice is to take No equal to the weighted geometric mean of the iV values occurring 
as data points.) As a consequence, the relative error on b is a factor of 3-4 smaller than in fits with 
No = 1. The error on A does not depend on the choice of N . 
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in which v has been fixed equal to 3/4 and the parameters a, b, A are free. In these fits we 
use only the data with A > some cutoff value A min ; we then vary A min systematically and 
investigate the quality of the fit (see Tables 13, 14 and 15). For A min > 60 the x 2 values 
are reasonable. The fits are stable and the error bars reasonable up to iV min ps 700; 
after this, the error bars increase drastically and the estimates should be considered 
unreliable. 

Let us consider first (R 2 ) and (R 2 g ), the radii for which we have the best statistics. 
The fit for the radius of gyration is extremely stable 9 and clearly suggests A = 1. 
No sub leading exponent with A < 1 appears to be present; in particular, this excludes 
A = 11/16 = 0.6875 unless the corresponding amplitude is extremely small. By contrast, 
the fit of {R 2 ) gives estimates that vary with A min : the estimate of A at first increases 
with N min , then flattens off at A 0.9. The theoretical prediction A = 11/16 is 
again excluded, but in this case a non-analytic correction A 1 < 1 is still possible a 
priori. However, we believe that a subleading exponent A 1 pa 0.9 is unlikely: after 
all, it does not agree with any theoretical prediction; and the observed behaviour can 
be explained equally well, as noted earlier, by the competition between two correction 
terms of opposite sign, provided that both terms are still sizable in the range of iV that 
we are considering. The fact that a range 500 < A < 4000 is insufficient to see clearly 
that A = 1 shows how difficult is the determination of A and explains the wide range 
of contradictory results found in previous work. 

We can also consider (R^)- The behaviour is similar to that observed for (R 2 ), 
although the errors are larger. Again, the data are compatible with A ps 0.9, but, as 
before, we believe that what we are observing is simply an effective exponent arising 
from the competition between two correction terms of opposite sign. Moreover, we are 
here probing shorter walks than in the case of (R 2 ), because of the lack of statistics for 
larger A: the fit is effectively dominated by the data in the range 300 ^ A < 1000. 

Since we have little evidence that Ai = 11/16, the most likely possibility is that 
A 1 = 3/2. We have therefore checked whether our data can be fitted by the simple 
Ansatz 

(R 2 ) N = N 2v {a + bN~ l + ciV- 3/2 ), (4.2) 

with a, b, c free parameters. As before, we perform several fits using only data with 
A > Amin, varying A m i n systematically. The results are reported in Tables 16, 17 and 
18. The fit quality is good even for the lowest value of A min : additional corrections 
do not play much role for A > 60. Notice that for (R 2 ) the constant c is very small, 
explaining why the previous fit gave A = 1 essentially without corrections. For {R 2 ) 
and (R^), c is instead sizable and of opposite sign with respect to b, giving in fit (4.1) 
an effective exponent A < 1. 

The leading amplitudes are extremely stable and, using the data with A min = 200, 

9 The x" 2 of the fit is somewhat large, and the corresponding confidence level too small. Since 
the confidence level actually gets worse as N m - ln grows, the cause does not seem to be corrections to 
scaling. The most likely interpretation is that our error bars are, for some unknown reason, somewhat 
underestimated for large values of N . (Perhaps we have underestimated the integrated autocorrelation 
time by failing to include a sufficient amount of the tail of the autocorrelation function.) 
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we can estimate 



a e 


= 0.77121 ± 0.00004 


(4.3) 


b e 


= 1.17 ±0.05 


(4.4) 




= -2.8 ±0.7, 


(4.5) 


a g 


= 0.108207 ±0.000007 


(4.6) 


b 9 


= 0.108 ±0.010 


(4.7) 


C 9 


= 0.0 ±0.2, 


(4.8) 




= 0.33903 ±0.00004 


(4.9) 


b m 


= 0.46 ±0.05 


(4.10) 


Cm 


= -1.1 ±0.7. 


(4.11) 



where the error bars are 68% confidence limits. We can compare these results with the 
series estimates (3.27)-(3.29), and note that they agree well within quoted errors with 
one exception: a g differs by three error bars from one set of the corresponding series 
estimate. We consider the stated errors of this series estimate to be anomalously low, 
compared to all other error estimates, and thus not to be taken literally. 

We have also considered the universal amplitude ratios A N and B N . The most no- 
table feature of the raw data for A N (Table 12) is its nonmonotonicity: at first A N 
decreases, reaching a minimum at N as 130; then it increases. This immediately sug- 
gests the presence of two correction-to-scaling terms of opposite sign, in agreement with 
the analysis presented above. We have therefore analysed the ratios A N and B N by 
performing a fit of the form 

Q N = a + bN~ l + cN~ 3/2 . (4.12) 

The results are reported in Tables 19 and 20. Again the quality of the fits is quite good, 
and we obtain the final estimates (again we use conservatively the data for A^ min = 200) 

A = 0.140310 ±0.000011 (4.13) 
B = 0.439614 ±0.000050, (4.14) 

where the error bars are 68% confidence limits. 10 

Next we analysed the CSCPS combination (2.11). Conformal-invariance theory [29, 
30] predicts that linijv-^oo -Fjv = 0, and we confirm this prediction numerically to very 
high precision: see Table 21. Therefore, fx is expected to scale as N 2v ~ 6 where 5 is some 
subleading exponent. Our data are consistent with 5 = 3/2, although not so accurate 
as to establish it unambiguously; in other words, /jy appears to approach a nonzero 

10 Cardy and Mussardo [39] have used the form-factor method, applied to the exact ^-matrix of the 
massive 0(n) model, to derive the estimates A w 0.126 and B w 0.420. This is impressive accuracy for 
a first-principles theoretical calculation; the approximately 5% error is about what one expects from 
the one-particle approximation used in this computation. 
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constant as iV — > oo. This means, as noted earlier, that the 1/N correction is absent 
within our errors. A fit to a constant gives the results reported in Table 22. The results 
show an initial downward trend with iV min and then increase again. A stable region may 
be identified for N min > 250. For iV min = 250 we have 

/ = 0.68 ±0.14. (4.15) 

This is in agreement with, but less precise than, the result / = 0.79 ± 0.03 obtained 
from series analysis. 

Finally, we have tried to see whether our estimated values of A and B are con- 
sistent with simple rational values that satisfy the CSCPS formula. If we require the 
denominators to be < 1000, then the only possible values anywhere near our estimates 
are 

23 40 

A = ~ 0.1402439, B = — « 0.4395604 (4.16) 

164 '91 V ; 

and 

91 95 

A = « 0.1404321, B = ^ 0.4398148. (4.17) 

648 '216 v ; 

Our data are, however, precise enough to clearly exclude both guesses. We therefore 
conjecture that A and B do not take simple rational values, even though one particular 
linear combination of them does. 



4.3 Comparison to series-extrapolation predictions 

We can also directly compare our raw Monte Carlo data to the extrapolation formulas 
(3.27)-(3.29). For this purpose we compute 

X 2 = Y, ( R2mc ~ R2scAcs)2 , (4.18) 

MCdata a MC 

where R 2 MC is the Monte Carlo estimate, o~mc the corresponding error, and R 2 eries the 
prediction of the extrapolations (3.27)-(3.29). For (R 2 ,) we find that (3.27) describes the 
numerical data rather well. Indeed, \ 2 = 29 for 19 data points. The small remaining 
discrepancy is mainly due to the error on the coefficients. Indeed, if we use for the 
ferromagnetic part the results of the Monte Carlo fit with N min = 40 reported in Table 
16 — these estimates are compatible with the exact-enumeration ones reported in (3.27)— 
the x 2 drops to 15. Analogous discussion applies to (R^)- If we use (3.29) we obtain 
X 2 = 55 with 15 data points. But again it is enough to replace the ferromagnetic 
coefficients obtained using exact enumeration with those obtained in the Monte Carlo 
fit, see Table 18 — they are fully compatible — to have x 2 = 19 with 15 points. The 
situation is worse for (R 2 )- Using all data we obtain x 2 — 80 with 18 data points. Such 
a result does not improve significantly if we change the coefficients in (3.28) within error 
bars. This is related to the fact we have already noticed that the leading amplitude for 
(R 2 )n reported in (3.28) significantly differs from the Monte Carlo estimate obtained 
for any value of iV min . 
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5 Conclusions and open questions 



In this study of the SAW correction-to-scaling exponents, we have seen a consistent 
picture emerging, given independent support by both Monte Carlo and series analysis. 
We have presented compelling evidence that the first non-analytic correction term in the 
generating function for both SAWs and SAPs, as well as in several Euclidean-invariant 
metric properties, is A 1 = 3/2, as predicted by Nienhuis some 20 years ago [4,5]. We 
find no evidence for the presence of an exponent Ai = 11/16 in SAWs and SAPs on 
the square and triangular lattices. Our analysis of the interplay between dominant and 
subdominant correction-to-scaling terms also enables us to explain quantitatively why 
many earlier analyses gave incorrect conclusions, predicting exponents Ai < 1. For 
certain observables, we find pairs of correction terms of opposite sign that conspire to 
give effective exponents that are smaller than both of the individual exponents. Thus, 
corrections behaving as a/N + b/N 3 ^ 2 with ab < were incorrectly identified with a 
single correction term c/N Al with Ai < 1. 

Monte Carlo and series analysis turn out to complement each other well. Series 
provide a basis for calculating the amplitudes of several subdominant asymptotic terms, 
while the Monte Carlo data frequently provide greater accuracy for the estimate of the 
leading amplitudes. 

We have also studied the asymptotic behaviour of several non-Euclidean-invariant 
quantities. Their leading behaviour is characterized by a new exponent A nr . We find 
compelling evidence that A nr = 2v on the square lattice and A nr = Au on the triangular 
lattice, confirming the conjecture of [43,44]. We also computed the leading correction- 
to-scaling exponent in these observables, finding A x 0.5. We are unaware of any 
theoretical prediction for this quantity. 

We have also determined the dominant and subdominant exponents characterizing 
the "antiferromagnetic singularity" of the square lattice. These exponent predictions 
are for the most part new. 

We also tested the CSCPS relation lini/v^oo-^jv — [cf. (2.10)], which arises from 
conformal field theory [29,30]. Both our Monte Carlo and series work are completely 
consistent with the CSCPS relation. Further, we find that the 1/N correction term in 
F N is absent, so that F N ~ const x iV~ Al with Ai = 3/2. The absence of this analytic 
correction-to-scaling term implies a new amplitude relation (2.44). 

Finally, we remark that our numerical estimates for the universal amplitude ratios 
A and B are now so precise as to allow us to rule out the possibility that these are 
rational numbers with small integer denominators. Some other universal amplitude 
ratios include powers of 7r [38], so it is possible that A and B are combinations of n and 
rational numbers; but there is no a priori reason why powers of n should enter into the 
amplitude ratios A and B. 

Acknowledgments 

The authors wish to thank Giovanni Ferraro for many helpful discussions and Raghu 
Varadhan for interesting suggestions. We also wish to thank Debra Bennett- Wood and 



35 



Bin Li for helping to generate some of the series and Monte Carlo data, respectively; 
Aleks Owczarek for helpful comments on the manuscript; and Peter Grassberger for 
providing high-precision Monte Carlo data that have been included in the analysis. 
Finally, we wish to thank an anonymous referee for important comments concerning the 
renormalization-group explanation of the cancellation of the analytic corrections to the 
CSCPS relation. 

This research was supported in part by the Consiglio Nazionale delle Ricerche (S.C. 
and A. P.), the Australian Research Council (A.J.G. and I. J.), the U.S. National Sci- 
ence Foundation grants DMS-8911273, DMS-9200719, PHY-9520978, PHY-9900769, 
PHY-0099393 and PHY-0424082 (A.D.S.), the U.S. Department of Energy contract 
DE-FG02-90ER40581 (A.D.S.), the NATO Collaborative Research grant CRG 910251 
(S.C. and A.D.S.), and a grant from the New York University Research Challenge Fund 
(A.D.S.). Acknowledgment is also made to the donors of the Petroleum Research Fund, 
administered by the American Chemical Society, for partial support of this research un- 
der grants 21091-AC7 and 25553-AC7B-C. A.N.R. and I.J. are pleased to acknowledge 
computational facilities provided by both VPAC and APAC. 

References 

[1] P.G. DeGennes, Scaling Concepts in Polymer Physics (Cornell Univ. Press, Ithaca, 
NY, 1979). 

[2] J. des Cloizeaux and G. Jannink, Polymers in Solution: Their Modelling and Struc- 
ture (Oxford University Press, New York, 1990). 

[3] L. Schafer, Excluded Volume Effects in Polymer Solutions (Springer- Verlag, Berlin- 
New York, 1999). 

[4] B. Nienhuis, Phys. Rev. Lett. 49, 1062 (1982). 

[5] B. Nienhuis, J. Stat. Phys. 34, 731 (1984). 

[6] H. Saleur, J. Phys. A: Math. Gen. 20, 455 (1987). 

[7] J. Adler, J. Phys. A: Math. Gen. 16, L515 (1983). 

[8] Z.V. Djordjevic, I. Majid, H.E. Stanley and R.J. dos Santos, J. Phys. A: Math. 
Gen. 16, L519 (1983). 

[9] I. Majid, Z.V. Djordjevic and H.E. Stanley, Phys. Rev. Lett. 51, 143 (1983). 

[10] V. Privman, Physica A 123, 428 (1984). 

[11] A.J. Guttmann, J. Phys. A: Math. Gen. 17, 455 (1984). 

[12] D.C. Rapaport, J. Phys. A: Math. Gen. 18, L201 (1985). 

[13] T. Ishinabe, Phys. Rev. B 37, 2376 (1988). 



36 



[14] T. Ishinabe, Phys. Rev. B 39, 9486 (1989). 

[15] D. MacDonald, D.L. Hunter, K. Kelly and N. Jan, J. Phys. A: Math. Gen. 25, 1429 

(1992) . 

[16] A.R. Conway, I.G. Enting and A.J. Guttmann, J. Phys. A: Math. Gen. 26, 1519 

(1993) , hep-lat/9211062. 

[17] A.R. Conway and A.J. Guttmann, Phys. Rev. Lett. 77, 5284 (1996). 

[18] S.R. Shannon, T.C. Choy and R.J. Fleming, Phys. Rev. B 53, 2175 (1996), cond- 
mat/9510162. 

[19] I. Jensen and A.J. Guttmann, J. Phys. A: Math. Gen. 32, 4867 (1999), cond- 
mat/9905291. 

[20] S. Havlin and D. Ben-Avraham, Phys. Rev. A 27, 2759 (1983). 
[21] J.W. Lyklema and K. Kremer, Phys. Rev. B 31, 3182 (1985). 

[22] D.L. Hunter, N. Jan and D. MacDonald, J. Phys. A: Math. Gen. 19, L543 (1986). 
[23] P.M. Lam, Phys. Rev. B 42, 4447 (1990). 

[24] I. Guim, H.W.J. Blote and T.W. Burkhardt, J. Phys. A: Math. Gen. 30, 413 (1997). 

[25] A.R. Conway and A.J. Guttmann, J. Phys. A: Math. Gen. 26, 1535 (1993), hep- 
lat/9211063. 

[26] I. Jensen, J. Phys. A: Math. Gen. 34, 7979 (2001). 

[27] S. Caracciolo, M.S. Causo, P. Grassberger and A. Pelissetto, J. Phys. A: Math. 
Gen. 32, 2931 (1999), cond-mat/9812267. 

[28] I. Jensen and A.J. Guttmann, J. Phys. A: Math. Gen. 33, L257 (2000). 

[29] J.L. Cardy and H. Saleur, J. Phys. A: Math. Gen. 22, L601 (1989). 

[30] S. Caracciolo, A. Pelissetto and A.D. Sokal, J. Phys. A: Math. Gen. 23, L969 (1990). 

[31] N. Madras and G. Slade, The Self-Avoiding Walk (Birkhauser, Boston-Basel- 
Berlin, 1993). 

[32] T. Hara and G. Slade, Commun. Math. Phys. 147, 101 (1992). 
[33] T. Hara and G. Slade, Reviews Math. Phys. 4, 235 (1992). 

[34] T. Hara, G. Slade and A.D. Sokal, J. Stat. Phys. 72, 479 (1993), hep-lat/9302003; 
78, 1187 (E) (1995). 

[35] A. Pelissetto and E. Vicari, Phys. Reports 368, 549 (2002), cond-mat/0012164. 



37 



[36] G.F. Lawler, O. Schramm and W. Werner, in Fractal Geometry and Applications: 
A Jubilee of Benoit Mandelbrot, Part 2, Proceedings of Symposia in Pure Mathe- 
matics #72 (American Mathematical Society, Providence RI, 2004), pp. 339-364, 
math.PR/0204277. 

[37] B. Li, N. Madras and A.D. Sokal, J. Stat. Phys. 80, 661 (1995), hep-lat/9409003. 

[38] J.L. Cardy and A.J. Guttmann, J. Phys. A: Math. Gen. 26, 2485 (1993), cond- 
mat/9303035. 

[39] J.L. Cardy and G. Mussardo, Nucl. Phys. B 410, 451 (1993), hep-th/9306028. 

[40] A.J. Guttmann and Y.S. Yang, J. Phys. A: Math. Gen. 23, L117 (1990). 

[41] A. L. Owczarek, T. Prellberg, D. Bennett- Wood and A.J. Guttmann, J. Phys. A: 
Math. Gen. 27, L919 (1994). 

[42] S. Caracciolo, M.S. Causo, and A. Pelissetto, J. Chem. Phys. 112, 7693 (2000), 
hep-lat/9910016. 

[43] M. Campostrini, A. Pelissetto, P. Rossi and E. Vicari, Phys. Rev. E 57, 184 (1998), 
cond-mat/9705086. 

[44] M. Campostrini, A. Pelissetto, P. Rossi, and E. Vicari, Phys. Rev. B 54, 7301 
(1996), hep-lat/9603002. 

[45] M. Caselle and M. Hasenbusch, Nucl. Phys. B 579, 667 (2000), hep-th/9911216. 

[46] M. Caselle, M. Hasenbusch, A. Pelissetto and E. Vicari, J. Phys. A: Math. Gen. 
35, 4861 (2002), cond-mat/0106372. 

[47] H. Cheng and T.T. Wu, Phys. Rev. 164, 719 (1967). 

[48] J. Stephenson, J. Math. Phys. 5, 1009 (1964); 11, 420 (1970). 

[49] F.J. Wegner, Phys. Rev. B 5, 4529 (1972). 

[50] V. Privman, P.C. Hohenberg and A. Aharony, in Phase Transitions and Critical 
Phenomena, Vol. 14, ed. C. Domb and J.L. Lebowitz (Academic Press, London- 
San Diego, 1991). 

[51] B.G. Nickel, Macromolecules 24, 1358 (1991). 

[52] A. Aharony and M.E. Fisher, Phys. Rev. B 27, 4394 (1983). 

[53] S. Gartenhaus and W.S. McCullough, Phys. Rev. B 35, 3299 (1987). 

[54] S. Gartenhaus and W.S. McCullough, Phys. Rev. B 38, 11688 (1988). 

[55] J. Salas and A.D. Sokal, Universal amplitude ratios in the critical two-dimensional 
Ising model on a torus, cond-mat/9904038vl, Section 3.1. 



38 



[5?; 

[58 
[59 

[60 
[61 
[62 

[63 
[64 

[65 
[66 
[67 
[68 
[69 
[70 

[71 
[72 
[73 
[74 



W.P. Orrick, B. Nickel, A.J. Guttmann and J.H.H. Perk, J. Stat. Phys. 102 , 795 

(2001) , cond-mat/0103074. 

M.F. Sykes and M.E Fisher, Phys. Rev. Lett. 1, 321 (1958). 
M.E. Fisher, Philos. Mag. 7, 1731 (1962). 

M.E. Fisher, in Statistical Physics, Weak Interactions, Field Theory (Lectures in 
Theoretical Physics, vol. 7C), edited by Wesley E. Brittin (University of Colorado 
Press, Boulder, 1965), pp. 1-159. 

S.S.C. Burnett and S. Gartenhaus, Phys. Rev. B 47, 7944 (1993). 

A.J. Guttmann and S.G. Whittington, J. Phys. A: Math. Gen. 11, 729 (1978). 

M. Campostrini, A. Pelissetto, P. Rossi and E. Vicari, Phys. Rev. E 65, 066127 

(2002) , cond-mat/0201180. 

A.J. Guttmann and J. Wang, J. Phys. A: Math. Gen. 24, 3107 (1991). 

D.N. Bennett- Wood, Numerical Studies of Self-Avoiding Walks, Ph.D. thesis, Uni- 
versity of Melbourne (1998). 

I. Jensen, J. Phys. A: Math. Gen. 37, 5503 (2004), cond-mat/0404728. 
I. Jensen, J. Phys. A: Math. Gen. 36, 5731 (2003), cond-mat/0301468. 
I.G. Enting and A.J. Guttmann, J. Phys. A: Math. Gen. 25, 2791 (1992). 
A.N. Rogers, Ph.D. thesis, University of Melbourne (2003). 
I. Jensen, J. Stat. Mech., P10008 (2004), cond-mat/0409039. 

A.J. Guttmann, in Phase Transitions and Critical Phenomena, vol. 13, C. Domb 
and J. Lebowitz, eds. (Academic Press, London, 1989). 

A.J. Guttmann, J. Phys. A: Math. Gen. 20, 1839 (1987). 

A. J. Guttmann, Discrete Math. 217, 167 (2000). 
M. Lai, Molec. Phys. 17, 57 (1969). 

B. MacDonald, N. Jan, D.L. Hunter and M.O. Steinitz, J. Phys. A: Math. Gen. 18, 
2627 (1985). 



[75] N. Madras and A.D. Sokal, J. Stat. Phys. 50, 109 (1988). 



39 



A 7 



\cn(RI)n 



\{N + l) 2 c N (R 2 g ) N \(N + l)c N (R 2 J N 





1 


4 


1 


1 


1 


2 


12 


8 


14 


11 


3 


36 


41 


116 


74 


4 


100 


176 


722 


390 


5 


284 


679 


3887 


1801 


6 


780 


2452 


18508 


7537 


7 


2172 


8447 


82160 


29684 


8 


5916 


28120 


340180 


110796 


9 


16268 


91147 


1351555 


399375 


10 


44100 


289324 


5136194 


1391809 


11 


120292 


902721 


18989580 


4741466 


12 


324932 


2777112 


68082102 


15783154 


13 


881500 


8441319 


239338055 


51704949 


14 


2374444 


25398500 


822629240 


166550157 


15 


6416596 


75744301 


2786064872 


530165200 


16 


17245332 


224156984 


9274487688 


1666083296 


17 


46466676 


658855781 


30521878637 


5188200085 


18 


124658732 


1924932324 


99086541810 


15993447527 


19 


335116620 


5593580859 


318742922236 


48946213794 


20 


897697164 


16175728584 


1014076260686 


148574713674 


21 


2408806028 


46572304083 


3202213457395 


448343690109 


22 


6444560484 


133556779740 


10019907413348 


1343838723819 


23 


17266613812 


381611332725 


31158181454688 


4008314601988 


24 


46146397316 


1086759598120 


96149048417484 


11888228627772 


25 


123481354908 


3085406711831 


295142819123871 


35114454662483 


26 


329712786220 


8735073410100 


900066956153270 


103219276329251 


27 


881317491628 


24665061125667 


2732505731274220 


302350533278086 


28 


2351378582244 


69477665745896 


8248829001382526 


881974673999634 


29 


6279396229332 


195265123427301 


24804499283684685 


2564984750250567 


30 


16741957935348 


547633156505396 


74221205928683512 


7432620540208579 


31 


44673816630956 


1532838884952299 


221333034486502552 


21480821587356344 


32 


119034997913020 


4282540754311160 


657175409080839632 


61884062343185928 


33 


317406598267076 


11944032183124129 


1945418024966109721 


177867529730724713 


34 


845279074648708 


33257656763184556 


5736868444918797822 


509789937527302533 


35 


2252534077759844 


92461749453584977 


16872515769277148908 


1458110752331771118 


36 


□yyo / 4U4yy ij44iz 


ZO668OO81589089 / 20 


4y4!>>lby48 12022 115 10 


A 1 eilArin tj r i c c; 1 torn 
41000001 /00IDDI60DO 


37 


15968852281708724 


711610318376609453 


144602103638139094373 


11846978451964635723 


38 


42486750758210044 


1970232464253179804 


421517897111033526836 


33661784893928852621 


39 


113101676587853932 


5448222121256407587 


1226106614925695416296 


95488104372916704748 


40 


300798249248474268 


15048127109659424048 


3556566126265101309980 


270319420601 198687276 


41 


800381032599158340 


41516822124396623905 


10296602201097594942225 


7641 12980773483151309 


42 


2127870238872271828 


114420546244580495788 


29734383555237595018738 


2155935506456910539997 


43 


5659667057165209612 


315023293116319316107 


85715948159939093661428 


6074721940189140309702 


44 


15041631638016155884 


866485020069260644664 


246526159900360427960866 


17087792773773847037390 


45 


39992704986620915140 


2381096560500892770793 


707895621958240811004041 


48007613778509804802835 


46 


106255762193816523332 


6537456672967882139948 


2028423176500049998353720 


134668241475007682635679 


47 


282417882500511560972 


17933790994378821974707 


5803786650111498416967400 


377337842568931343367952 


48 


750139547395987948108 


49156721156019804756024 


16573780740476629162797384 


1055798043293664579046904 


49 


1993185460468062845836 


134634067180094086612595 


47265785843405796903827851 


2951082969829519054855011 


50 


5292794668724837206644 


368470259697878846546356 


134553477427137032959671590 


8237877051162727078333841 


51 


14059415980606050644844 


1007714451664000011164731 


382563667789601278749936172 


22973891188237620817265778 


52 


37325046962536847970116 


2754055291248183384820416 


1085908845295037795117686914 


63992592137176740914658486 


53 


99121668912462180162908 


7521747594983831934039415 


3078819835262093742193456615 


178090816717561120403807465 


54 


263090298246050489804708 


20529906173170669487082516 


8715811544642479874846017668 


495070634995239123112723883 


55 


698501700277581954674604 


55999859678542460919790667 


24647267299559828818119548152 


1375109866469718419312054508 


56 


1853589151789474253830500 


152661200274313840382160720 


69599920787167645242434336340 


3815534534588842152972448564 


57 


4920146075313000860596140 


415930558458671192311805699 


196344377133567516049914306371 


10578959386509767813473708087 


58 


13053884641516572778155044 


1132589111567068180634238436 


553157627020295519719444493138 


29302775940073964314770902241 


59 


34642792634590824499672196 


3082415161154176613766926049 


1556962657584619410795878310156 


81108542428478523911668978514 



Table 1: Exact enumeration data for SAWs on the square lattice. 
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Table 2: Exact enumeration data for SAWs on the triangular lattice. 



41 



N 



N 



jC N (r 6 ) N 



jC N (r 8 ) N 



1 


1 


1 


1 


2 


24 


80 


288 


3 


233 


1481 


10313 


4 


1552 


15584 


171712 


5 


8261 


118741 


1876421 


6 


40128 


761824 


15997248 


7 


174687 


4216895 


113009823 


8 


711744 


21139264 


699292800 


9 


2756691 


98246971 


3911019843 


10 


10258032 


430155712 


20197992960 


11 


36953225 


1794576465 


97801373081 


12 


129595552 


7194227712 


449049597184 


13 


444358551 


27891276903 


1971835847895 


14 


1494601312 


105092615072 


8336039677888 


15 


4944384005 


386372087101 


34107295967573 


16 


16121969312 


1390424839040 


135635987698688 


17 


51903980173 


4910490964373 


526075527334141 


18 


165229382704 


17055786755328 


1995781800318592 


19 


520720306083 


58367380590987 


7423576318235379 


20 


1626289219696 


197097871552608 


27128806075092160 


21 


5037880731363 


657614956490835 


97570424122840995 


22 


15491105783776 


2170327643009376 


345877992391828288 


23 


47313566966717 


7091919679833573 


1210056084152236397 


24 


143616941038800 


22964364302956192 


4182680502669028416 


25 


433471181567175 


73742760159367607 


14298788347356195303 


26 


1301492251611088 


234986541658461504 


48385725874370354944 


27 


3888842767461723 


743493757302422163 


162197418768856363467 


28 


11567743361677920 


2336936884325400320 


538986352207098913536 


29 


34265929488742837 


7300504880342236965 


1776583932652252100533 


30 


101107717070386128 


22676475755899170368 


5811758654106692557056 


31 


297251719690114411 


70061223151862034731 


18878148082088838145579 


32 


870928677714199072 


215380608263460514688 


60916738136365328424448 



Table 3: Exact enumeration data for SAWs on the square lattice. Rotationally invariant 
moments. 
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4284776312224 


1004647110417216 


21 


51978553251 


13337579454483 


3367181961982563 


22 


149793700032 


41203857605920 


11161461862577856 


23 


430013901309 


126414733312805 


36627023206822344 


24 


1230085625008 


385396000445280 


119090886021960640 


25 


3507275950151 


1168118922135351 


383951817561778304 


26 


9970080369360 


3521536905190720 


1228239234255697152 


27 


28262765992155 


10563701297658387 


3900781003385081163 


28 


79911109071584 


31542297194620416 


12305781139629052160 


29 


225398486017269 


93778343662150501 


38579685692175173877 


30 


634334283147728 


277695329500224576 


120249249844158963968 


31 


1781434024153067 


819226661598869419 


372773226138821611691 


32 


4993035148467488 


2408303159048790400 


1149726544611189212672 



Table 4: Exact enumeration data for SAWs on the square lattice. Non-rotationally- 
invariant moments. 



43 



N 


lc N (r 4 ) N 


lc N (r 6 ) N 


lc N (r 8 ) N 


1 


1 


1 


1 


9 




190 


49D 


o 


59Q 


3337 




4 


5454 


53358 




5 


461 fiQ 


fi^SI 85 


Q79584Q 


6 


344428 


6221884 


126456796 


7 


2352769 


53647549 


1380311377 


8 


15060090 


420194610 


13256099610 


9 


91701871 


3057404227 


115436446639 


10 


536695548 


20985811596 


930424151244 


11 


3041620465 


137386509145 


7043622904369 


12 


16784388968 


864860517248 


50628426215432 


13 


90564392107 


5267777416675 


348381660817291 


14 


479388030946 


31195330864090 


2309752803978322 


15 


2495911050345 


180302936975925 


14829889785198921 


16 


12808562012852 


1020273312831596 


92588968138883348 


17 


64901247920059 


5666785971562159 


564018088195524619 


18 


325170745810666 


30958613694252346 


3361656410695492858 


19 


1612871078099977 


166655904211475269 


19649810153900928217 


20 


7927964773508104 


885333404431705216 


112867935343523712424 


21 


38652473796950531 


4647266177110051355 


638156051853903239891 


22 


187056802703356296 


24130876287242419704 


3556820174217345377400 



Table 5: Exact enumeration data for SAWs on the triangular lattice. Rotationally 
invariant moments. 



44 



TV 


g l-jY \ ' w o ui/ j j\j 


-Cat(t 8 COS 66) at 
g '-'iV \ ^-wo \jyj j j\j 


1 


1 


1 


2 


12 


96 


3 


97 


1609 


4 


654 


17454 


5 


3977 


151649 


6 


22684 


1149148 


7 


123721 


7935829 


8 


652842 


51236610 


9 


3357439 


314319571 


10 


16914348 


1852261068 


11 


83777857 


10566164665 


12 


409089560 


58677117008 


13 


1973505067 


318573263587 


14 


9421326322 


1696583222746 


15 


44567944521 


8885936605365 


16 


209144745044 


45868003547852 


17 


974497840243 


233746526628199 


18 


4511869867210 


1177691322037546 


19 


20770914530257 


5873419831448317 


20 


95130303643048 


29024633960838784 


21 


433664585252891 


142245181072370291 


22 


1968525488778840 


691879404495232056 



Table 6: Exact enumeration data for SAWs on the triangular lattice. Non-rotationally- 
invariant moments. 
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19 





712401 





847374 


-1 


932382 


1 


872489 


20 





712338 





854132 


-1 


970591 


1 


933229 


21 





712262 





862841 


-2 


021181 


2 


015860 


22 





712207 





869360 


-2 


060036 


2 


080978 


23 





712154 





876015 


-2 


100682 


2 


150783 


24 





712108 





882131 


-2 


138917 


2 


217997 


25 





712067 





887809 


-2 


175211 


2 


283237 


26 





712030 





893200 


-2 


210408 


2 


347867 


27 





711996 





898252 


-2 


244079 


2 


410972 


28 





711965 





903032 


-2 


276568 


2 


473071 


29 





711937 





907558 


-2 


307916 


2 


534132 


30 





711912 





911851 


-2 


338197 


2 


594197 


31 





711888 





915933 


-2 


367500 


2 


653354 


32 





711867 





919820 


-2 


395884 


2 


711636 


33 





711847 





923529 


-2 


423411 


2 


769094 


34 





711828 





927073 


-2 


450135 


2 


825769 


35 





711811 





930464 


-2 


476105 


2 


881700 


36 





711795 





933713 


-2 


501366 


2 


936922 


37 





711780 





936831 


-2 


525957 


2 


991468 


38 





711766 





939827 


-2 


549918 


3 


045367 


39 





711753 





942708 


— 2 


573280 


3 


098647 


40 





711741 





945482 


-2 


596077 


3 


151335 


19 





099915 





050034 





335467 


-0 


323304 


20 





099919 





049550 





338206 


-0 


327657 


21 





099919 





049529 





338324 


-0 


327850 


22 





099920 





049454 





338770 


-0 


328598 


23 





099919 





049539 





338252 


-0 


327708 


24 





099918 





049696 





337275 


-0 


325990 


25 





099917 





049907 





335921 


-0 


323556 


26 





099915 





050176 





334165 


-0 


320333 


27 





099913 





050482 





332124 


-0 


316508 


28 





099911 





050822 





329820 


-0 


312103 


29 





099908 





051186 





327296 


-0 


307186 


30 





099906 





051570 





324587 


-0 


301813 


31 





099904 





051969 





321720 


-0 


296026 


32 





099901 





052380 





318722 


-0 


289869 


33 





099899 





052799 





315611 


-0 


283376 


34 





099897 





053224 





312408 


-0 


276582 


35 





099895 





053652 





309126 


-0 


269516 


36 





099893 





054083 





305781 


-0 


262203 


37 





099891 





054513 





302384 


-0 


254667 


38 





099889 





054943 





298944 


-0 


246930 


39 





099887 





055372 





295472 


— 


239010 


40 





099885 





055797 





291974 


-0 


230925 


19 





313864 





190168 


-0 


037972 






20 





313817 





192786 


-0 


045478 






21 





313772 





195399 


-0 


053166 






22 





313731 





197964 


-0 


060906 






23 





313691 





200488 


-0 


068707 






24 





313655 





202966 


-0 


076542 






25 





313620 





205394 


-0 


084387 






26 





313588 





207772 


-0 


092232 






27 





313558 





210098 


-0 


100063 






28 





313529 





212373 


-0 


107868 






29 





313502 





214597 


-0 


115640 






30 





313477 





216769 


-0 


123371 






31 





313453 





218892 


-0 


131055 






32 





313430 





220965 


-0 


138688 






33 





313409 





222990 


-0 


146266 






34 





313388 





224969 


-0 


153785 






35 





313369 





226902 


-0 


161244 






36 





313351 





228791 


-0 


168639 






37 





313333 





230637 


-0 


175971 






38 





313317 





232441 


-0 


183237 






39 





313301 





234205 


— 


190437 






40 





313286 





235930 


— 


197570 







Table 7: Fit to (R 2 ) = N 3 ^ 2 (a + ^N' 1 + a 2 N~ 3 ^ + a 3 N~ 2 ) for SAWs on the triangular 
lattice. Data is for (R 2 ) N at top, then (R 2 ) N at middle, then (R^) n at bottom. The fit 
for (R 2 m )N includes only terms up to order (N~ 3 / 2 ). 
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N 




ao,e 




Ol,e 




d2,e 






















33 


0. 


.771324 


1. 


.127229 


-2. 


,737952 


3. 


.927648 


-2. 


.772841 


0. 


.124654 


-0 


,037820 





,655979 


34 


0. 


.771291 


1. 


.137155 


-2. 


.846606 


4. 


373539 


-3 


.423234 


0. 


.124636 


-0 


.036728 





.639632 


35 


0. 


.771305 


1. 


.132833 


-2 


,798505 


4. 


.172870 


-3 


.125666 


0. 


.124628 


-0 


,036255 





,632312 


36 


0. 


.771288 


1, 


.138263 


-2. 


.859903 


4. 


433123 


-3 


.517795 


0. 


.124619 


-0 


,035663 





.622865 


37 


0. 


.771293 


1. 


136612 


-2. 


.840943 


4. 


.351507 


-3. 


.392907 


0. 


.124616 


-0. 


,035484 





.619916 


38 


0. 


.771279 


1, 


.141331 


-2. 


.895949 


4. 


.591870 


-3. 


766261 


0. 


.124609 


-0 


,034975 





.611270 


39 


0. 


.771286 


1. 


.138877 


-2. 


.866927 


4. 


.463195 


-3. 


.563461 


0. 


.124605 


-0 


,034711 





,606663 


40 


0. 


.771271 


1. 


.144437 


-2. 


933613 


4. 


763067 


-4. 


.042810 


0. 


.124596 


-0 


,034116 





.595971 


41 


0. 


.771280 


1. 


.141243 


-2. 


,894778 


4. 


586020 


-3. 


.755880 


0. 


.124592 


-0 


.033776 





.589685 


42 


0. 


.771265 


1. 


.146886 


-2. 


,964317 


4. 


.907314 


-4. 


.283596 


0. 


.124584 


-0 


.033177 





,578325 


43 


0. 


.771273 


1. 


.143785 


-2. 


925602 


4. 


726095 


-3. 


.982046 


0. 


.124580 


-0 


.032849 





.571943 


44 


0. 


.771260 


1. 


.149042 


-2. 


.992074 


5. 


.041219 


-4. 


.513116 


0. 


.124573 


-0 


.032296 





.560888 


45 


0. 


.771267 


1. 


.146099 


-2. 


954406 


4. 


.860421 


-4. 


,204625 


0. 


.124569 


-0 


.031987 





.554570 


46 


0. 


.771255 


1. 


.151169 


-3 


.020088 


5. 


.179507 


-4. 


,755692 


0. 


.124563 


-0 


,031458 


0. 


.543460 


47 


0. 


.771262 


1. 


.148199 


-2. 


.981150 


4. 


.988102 


-4. 


,421211 


0. 


.124559 


-0 


,031148 


0. 


.536820 


48 


0. 


.771250 


1, 


.153164 


-3 


.047007 


5. 


.315582 


-5. 


,000119 


0. 


.124553 


-0 


,030633 





.525499 


49 


0. 


.771257 


1. 


.150204 


-3 


.007298 


5. 


.115889 


-4. 


,643109 


0. 


.124550 


-0 


.030327 


0. 


.518620 


50 


0. 


.771246 


1. 


.155004 


-3 


.072388 


5. 


.446844 


-5. 


.241340 


0. 


.124545 


-0. 


,029833 


0. 


.507257 


51 


0. 


.771253 


1. 


.152113 


-3. 


,032759 


5. 


.243165 


-4. 


.869182 


0. 


.124541 


-0 


,029536 





,500287 


52 


0. 


.771243 


1. 


.156733 


-3. 


,096758 


5, 


.575580 


-5. 


.483004 


0. 


.124536 


-0 


.029063 





.488950 


53 


0. 


.771249 


1. 


.153911 


-3. 


,057261 


5, 


368299 


-5. 


096275 


0. 


.124534 


-0 


.028775 


0. 


.481902 


54 


0. 


.771240 


1. 


.158369 


-3. 


120296 


5. 


.702470 


-5. 


726090 


0. 


.124529 


-0 


.028322 





.470577 


55 


0. 


.771245 


1. 


.155612 


-3. 


,080928 


5, 


.491689 


-5. 


.324869 


0. 


.124526 


-0 


.028042 


0. 


.463455 


56 


0. 


.771237 


1. 


.159915 


-3. 


.142971 


5. 


.827124 


-5. 


969610 


0. 


.124522 


-0 


.027607 





,452155 


57 


0. 


.771242 


1. 


.157229 


-3. 


.103869 


5. 


.613694 


-5. 


.555442 


0. 


.124520 


-0 


.027337 





,444986 


58 


0. 


.771234 


1. 


.161379 


-3. 


.164847 


5. 


.949657 


-6. 


.213520 


0. 


.124516 


-0 


.026920 





,433734 


59 


0. 


.771239 


1. 


.158766 


-3. 


.126105 


5. 


734233 


-5. 


.787657 


0. 


124513 


-0 


.026658 





.426540 



Table 8: Fit (R 2 e ) N = iV 3 / 2 (a , e + a^iV" 1 + a 2:e N~ 3 / 2 + a^ e N~ 2 + a 4 , e N' 5 / 2 ) + 
(_l)A f iV- 11 /32( 6o e + b 1;e /N + b 2 , e /N2) for SAWs on the square lattice. 
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N 




°0,9 




°1,<? 




a 2,g 




Kg 




Kg 


33 


0. 


108213 


0. 


103702 


0. 


100792 





,008333 


0. 


,004174 


34 


0. 


108216 


0. 


103450 


0. 


101736 





,008370 


0. 


,003016 


35 


0. 


108218 


0. 


103242 


0. 


102526 





,008339 


0. 


,003994 


36 


0. 


108220 


0. 


103040 


0, 


103306 





,008369 





,003017 


37 


0. 


108221 


0. 


102873 


0. 


103959 


0. 


.008345 


0. 


.003842 


38 


0. 


108223 


0. 


102712 


0, 


104599 


0. 


.008368 


0. 


.003025 


39 


0. 


108224 


0, 


102579 


0, 


105134 


0. 


.008349 


0. 


.003715 


40 


0. 


108225 


0. 


102452 


0, 


105653 


0. 


.008367 


0. 


.003039 


41 


0. 


108226 


0, 


102347 


0. 


106085 


0. 


.008352 


0. 


.003607 


42 


0, 


108227 


0. 


102249 


0, 


106500 


0. 


.008366 


0. 


,003058 


43 


0. 


108228 


0, 


102168 


0. 


106842 


0. 


.008355 





,003515 


44 


0. 


108228 


0. 


102092 


0. 


107166 


0. 


.008365 


0. 


,003078 


45 


0. 


108229 


0, 


102032 


0, 


107429 





.008357 





,003436 


zLfi 


u. 


luozzy 


n 
u. 


1 m Q7R 
iu iy / u 


u. 




n 
u 


.UUoOUO 


n 
u 


UUO 1UU 


47 


0. 


108230 


0. 


101933 


0, 


107867 


0. 


.008359 





,003367 


48 


0. 


108230 


0. 


101894 


0. 


108043 





.008364 


0. 


.003122 


49 


0. 


108230 


0. 


101865 


0. 


108175 





.008360 


0. 


.003307 


50 


0. 


108230 


0. 


101840 


0. 


108290 


0. 


.008363 


0. 


.003145 


51 


0. 


108230 


0. 


101823 


0. 


108367 





.008361 


0. 


.003254 


52 


0. 


108230 


0. 


101810 


0. 


108428 





.008363 


0. 


.003166 


53 


0. 


108230 


0. 


101804 


0. 


108456 





.008362 


0. 


.003207 


54 


0. 


108230 


0. 


101801 


0. 


108470 





,008363 


0. 


,003187 


55 


0. 


108230 


0. 


101804 


0, 


108455 





,008363 


0. 


,003165 


56 


0. 


108230 


0. 


101810 


0. 


108427 





,008362 


0. 


,003207 


57 


0. 


108230 


0. 


101821 


0. 


108373 





,008364 


0. 


,003128 


58 


0. 


108230 


0. 


101834 


0. 


108307 





,008362 





,003226 


59 


0. 


108230 


0. 


101852 


0. 


108219 





,008364 





,003094 



Table 9: Fit (R 2 g ) N = iV 3 / 2 (a , 9 + a^iV" 1 + a^iV" 3 / 2 ) + (-1) N N' 11 / 32 ^ + b hg /N) for 
SAWs on the square lattice. 
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N 








0>l,m 




0>2,m 














33 





.339223 





.401114 


-0 


.849196 


1 


.011789 


0. 


.036366 


0. 


,020608 


34 





.339252 





.395837 


-0 


.810037 





.930082 


0. 


.036379 


0. 


.020222 


35 


0. 


.339210 





.403819 


-0 


.870217 


1 


.057665 


0. 


.036397 


0. 


,019642 


36 





.339233 





.399343 


-0 


.835945 





.983880 


0. 


.036407 


0. 


,019318 


37 


0. 


.339198 





.406364 


-0 


.890517 


1 


.103143 


0. 


.036422 


0. 


,018812 


38 





.339216 





.402541 


-0 


.860364 


1 


.036281 


0. 


.036430 





,018538 


39 





.339188 





.408723 


-0 


.909807 


1 


.147474 


0. 


.036442 


0. 


,018096 


40 





.339202 





.405479 


-0 


.883503 


1 


.087501 


0. 


.036448 


0. 


,017866 


41 





.339178 





.410928 


-0 


.928281 


1 


.190967 





.036458 


0. 


.017480 


42 





,339190 





.408177 


-0 


.905379 


1 


.137357 





.036463 





.017286 


43 





.339170 





.413005 


-0 


.946085 


1 


.233861 


0. 


.036472 


0. 


.016947 


44 





.339180 





.410665 


-0 


.926109 


1 


.185914 





.036476 





.016783 


45 





.339163 





.414963 


-0 


.963242 


1 


.276125 


0. 


.036483 


0. 


.016483 


46 


0. 


.339170 





.412973 


-0 


.945843 


1 


.233355 





,036486 





,016345 


47 


0. 


.339156 





.416811 


-0 


.979781 


1 


.317746 





,036492 





,016080 


48 


0. 


.339162 





.415122 


-0 


.964678 


1 


.279766 





,036495 





,015963 


49 


0. 


.339150 





.418559 


-0 


.995753 


1 


.358778 





,036500 





,015727 


50 





.339155 





.417130 


-0 


.982689 


1 


.325201 





,036502 





,015629 


51 





.339144 





.420217 


-1 


.011200 


1 


.399256 


0. 


.036506 


0. 


.015418 


52 





,339148 





.419011 


-0 


.999949 


1 


.369730 


0. 


.036508 


0. 


.015336 


53 


0. 


.339139 





.421791 


-1 


.026155 


1 


.439200 


0. 


.036512 


0. 


.015148 


54 





.339142 





.420779 


-1 


.016525 


1 


.413417 


0. 


.036513 


0. 


.015079 


55 





.339134 





.423289 


-1 


.040647 


1 


.478633 


0. 


.036516 


0. 


.014910 


56 


0. 


.339137 





.422446 


-1 


.032471 


1 


.456316 


0. 


.036517 


0. 


,014854 


57 





.339130 





.424716 


-1 


.054707 


1 


.517578 


0. 


.036520 


0. 


,014702 


58 


0. 


.339132 





.424021 


-1 


.047836 


1 


.498474 


0. 


.036521 


0. 


.014655 


59 





.339126 





.426078 


-1 


.068359 


1 


.556056 


0. 


.036523 





.014518 



Table 10: Fit (R 2 J N = iV 3 / 2 (a , m + a^iV" 1 + a 2 ^N~^ 2 + a^ m N~ 2 ) + 
(-l) Ar A^- 11 /32(ft 0m + b hm /N) for SAWs on the square lattice. 
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N 








TO. 


V 








Fn{RI)n 


40 


200.106 


± 


0. 


.011 
























fin 


364.977 


± 


0. 


.012 


51.2026 


± 


U 


uuo i 


160.10 


± 


0, 


,04 


0.71 


± 


0.08 


80 


559.656 


± 


0. 


.033 


78.4660 


± 


0. 


.0051 


245.51 


± 


0, 


06 


0.91 


± 


0.13 


100 


780.245 


± 


0. 


.028 


109.351 


± 


0. 


.007 


342.34 


± 


0, 


09 


1.06 


± 


0.18 


120 


1023.812 


± 


0. 


.064 


143.488 


± 


0. 


.010 


449.39 


± 


0. 


,12 


1.02 


± 


0.24 


140 


1288.631 


± 


0. 


.064 


180.599 


± 


0. 


.013 


565.93 


± 


0, 


.16 


0.67 


± 


0.32 


150 


1428.367 


± 


0. 


.035 


200.185 


± 


0. 


.015 


627.51 


± 


0. 


.18 


0.32 


± 


0.36 


180 


1875.399 


± 


0. 


.081 


262.858 


± 


0. 


.020 


823.95 


± 


0, 


,13 


0.38 


± 


0.26 


200 


2195.00 


± 


0. 


.11 


307.676 


± 


0. 


.024 


964.36 


± 


0, 


,13 


0.52 


± 


0.27 


250 


3064.11 


± 


0. 


.11 


429.499 


± 


0. 


.024 


1346.26 


± 


0, 


,19 


0.59 


± 


0.39 


300 


4024.66 


± 


0. 


.17 


564.226 


± 





.021 


1768.44 


± 


0, 


,07 


0.71 


± 


0.18 


400 


6190.01 


± 





.26 


867.890 


± 


0. 


.033 


2720.27 


± 


0, 


,12 


0.62 


± 


0.28 


500 


8645.61 


± 





.31 


1212.26 


± 





.07 


3799.66 


± 


0, 


,25 


0.57 


± 


0.57 


700 


14311.70 


± 





.57 


2007.10 


± 





.11 
















1000 


24421.12 


± 


0. 


.85 


3425.28 


± 


0. 


.11 


10734.5 


± 


0, 


53 


1.15 


± 


1.20 


1400 


40439.29 


± 


2 


.11 


5671.95 


± 





,42 
















2000 


69028.33 


± 


3 


.82 


9684.17 


± 





,79 


30347.0 


± 


3, 


,5 


-0.6 


± 


7.5 


3000 


126789.4 


± 


7 


.9 


17790.3 


± 


2 


,1 
















4000 


195162.3 


± 


12.2 


27376.1 


± 


2 


,7 


85773.6 


± 


13. 


40 


± 


27 



Table 11: Our Monte Carlo data for radii (R 2 )n as a function of walk length N. Errors 
(one standard deviation) are shown in parentheses. 
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A T 




A N 






7) 






Oat 






F N 


60 





140290 


(10) 





A O O f A A 

438644 


(100) 





319826 


(76) 





r\r\ i f \ r" n i c\r\r\\ 

001957 (200) 


80 





1 A r\c\r\ A 

140204 


(12) 





A O O Z" 1 O *"l 

438689 


(120) 





319598 


(85) 





001635 (230) 


100 





140149 


(11) 





438756 


(120) 





319425 


(88) 





001354 (240) 


120 





140150 


(13) 





A O O /"I O O 

438938 


(120) 





319294 


/aa\ 

(88) 





000992 (240) 


140 





-1 A r\ -1 A o 

140148 


(12) 





439171 


(130) 





319118 


(93) 





a a a r 1 r\ f r\v r\\ 

000518 (250) 


150 





140149 


(11) 





439320 


(130) 





319014 


(94) 





000225 (250) 


180 





140161 


(12) 





439347 


(71) 





319022 


(55) 





r\r\r\c\r\ A / -i A r\ \ 

000204 (140) 


200 





140171 


(13) 





439344 


(63) 





319047 


/ A r\\ 

(49) 





000236 (120) 


250 





140171 


(9) 





A O AO r 

439365 


(64) 





319030 


(48) 





000193 (130) 


300 





140192 


(8) 





439402 


(26) 





319052 


/-, n \ 

(18) 





000177 (45) 




n 
u 




K°) 


n 






n 
u 






n 
u 


UUU1U1 V / 


500 





140216 


(10) 





439490 


(33) 





319043 


(29) 





000066 (66) 


700 





140242 


(10) 


















1000 





140259 


(7) 





439557 


(27) 





319091 


(19) 





000047 (48) 


1400 





140258 


(13) 


















2000 





140293 


(14) 





439630 


(56) 





319115 


(45) 




0.000008 (110) 


3000 





140314 


(19) 


















4000 





140274 


(16) 





439499 


(71) 





319167 


(57) 





000204 (140) 



Table 12: Our Monte Carlo data for amplitude ratios as a function of walk length. Errors 
(one standard deviation) are shown in parentheses. 



N ■ 

1 v mm 


a e 


b e 


A 


x 2 


DF 


CL 


40 


0.770998(23) 


0.001644(16) 


0.855(3) 


40.9 


16 


0.1% 


60 


0.771054(25) 


0.001589(19) 


0.869(4) 


14.9 


15 


46.2% 


80 


0.771081(30) 


0.001559(26) 


0.878(7) 


12.2 


14 


59.0% 


100 


0.771106(32) 


0.001532(29) 


0.887(8) 


8.4 


13 


82.0% 


120 


0.771122(37) 


0.001512(36) 


0.894(12) 


7.6 


12 


81.8% 


140 


0.771140(38) 


0.001491(38) 


0.902(13) 


4.8 


11 


93.9% 


150 


0.771140(39) 


0.001492(39) 


0.902(13) 


4.8 


10 


90.2% 


180 


0.771160(45) 


0.001467(47) 


0.914(19) 


4.1 


9 


90.7% 


200 


0.771148(49) 


0.001481(54) 


0.906(23) 


3.7 


8 


88.1% 


250 


0.771155(55) 


0.001472(62) 


0.911(29) 


3.6 


7 


82.0% 


300 


0.771148(64) 


0.001481(73) 


0.905(39) 


3.6 


6 


73.2% 


400 


0.771136(82) 


0.001495(94) 


0.895(58) 


3.5 


5 


61.8% 


500 


0.771161(93) 


0.001469(104) 


0.921(80) 


3.3 


4 


51.0% 


700 


0.771249(110) 


0.001396(109) 


1.051(156) 


2.3 


3 


51.1% 


1000 


0.770975(410) 


0.001571(344) 


0.692(332) 


0.9 


2 


62.4% 



Table 13: Fit (R 2 e ) = N 2v [a e + & e (iV/750) A ]. DF is the number of degrees of freedom 
and CL is the confidence level of the fit. 
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N ■ 

1 ' mm 


a g 


h 


A 


V 2 
A 


DF 

1 J X 


C!L 


fin 

w w 


108213(4') 


000142(3") 


1 037(8') 


19 9 

X .J . .J 


15 


17 5% 

X 1 . W /(J 


sn 

OW 


W. XWOZjWC/IW 1 


D 0001 46(4") 

W.WWWXtzUIt: 1 


1 093(1 1 ") 

X . WZjO 111 1 


1 fi 8 


14 

11 


2fi 7% 


1 nn 

X WW 


108205(5^1 


0001 51 (5") 

W.WWWXWXIW I 


1 007(1 5") 

X . WW 1 1 lu 1 


14 2 


1 3 

X(J 


35 9% 

Ua . <v /(J 


1 2D 


108205(6^ 


0001 50(5^ 

W.WWWXWW1W ) 


1 008(18^ 

X . UU(J 1 X (J J 


14 2 


1 2 

X Zj 


28 8% 

ZjO .O /(J 


140 


0.108207(6) 


0.000148(6) 


1.018(22) 


13.5 


11 


25.9% 


150 


0.108207(7) 


0.000148(7) 


1.019(26) 


13.5 


10 


19.5% 


180 


0.108210(8) 


0.000145(8) 


1.033(34) 


13.1 


9 


15.7% 


200 


0.108208(9) 


0.000147(9) 


1.023(41) 


12.9 


8 


11.4% 


250 


0.108205(10) 


0.000151(11) 


1.000(53) 


12.5 


7 


8.5% 


300 


0.108208(11) 


0.000147(13) 


1.025(66) 


12.1 


6 


6.0% 


400 


0.108201(15) 


0.000155(17) 


0.966(100) 


11.5 


5 


4.2% 


500 


0.108198(19) 


0.000158(21) 


0.937(146) 


11.5 


4 


2.2% 


700 


0.108213(20) 


0.000146(19) 


1.161(286) 


10.6 


3 


1.4% 


1000 


0.107954(1213) 


0.000384(1200) 


0.206(779) 


8.8 


2 


1.3% 



Table 14: Fit (R 2 g ) = N 2u [a g + b g (N/750) A ]. DF is the number of degrees of freedom 
and CL is the confidence level of the fit. 



N ■ 

1 v mm 






A 


x 2 


DF 


CL 


60 


0.338957(32) 


0.000635(27) 


0.850(17) 


15.4 


12 


21.9% 


80 


0.338932(38) 


0.000661(33) 


0.828(23) 


13.3 


11 


27.1% 


100 


0.338940(41) 


0.000652(38) 


0.836(29) 


13.1 


10 


21.7% 


120 


0.338971(43) 


0.000618(41) 


0.871(36) 


10.2 


9 


33.6% 


140 


0.339008(44) 


0.000577(42) 


0.917(42) 


5.7 


8 


68.1% 


150 


0.339018(45) 


0.000566(44) 


0.930(47) 


5.3 


7 


62.7% 


180 


0.339013(48) 


0.000571(48) 


0.923(52) 


5.2 


6 


52.1% 


200 


0.338990(56) 


0.000598(58) 


0.890(62) 


4.2 


5 


51.7% 


250 


0.338962(71) 


0.000630(75) 


0.852(79) 


3.7 


4 


45.3% 


300 


0.338955(75) 


0.000638(81) 


0.842(84) 


3.6 


3 


31.3% 


400 


0.338915(115) 


0.000683(126) 


0.784(136) 


3.3 


2 


19.2% 


500 


0.338896(146) 


0.000702(155) 


0.752(186) 


3.2 


1 


7.2% 



Table 15: Fit (R 2 m ) = N 2u [a m + b m (N/750) A ]. DF is the number of degrees of freedom 
and CL is the confidence level of the fit. 
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-^min 


a e 


b e 


c e 


x 2 


DF 


CL 


40 


0.771261(16) 


1.093(7) 


-1.93(4) 


15.4 


16 


49 


4% 


60 


0.771247(18) 


1.104(9) 


-2.03(7) 


12.1 


15 


67 


1% 


80 


0.771224(20) 


1.131(15) 


-2.30(14) 


6.7 


14 


94 


6% 


100 


0.771224(22) 


1.131(18) 


-2.30(18) 


6.7 


13 


91 


7% 


120 


0.771211(25) 


1.150(26) 


-2.52(29) 


5.7 


12 


93 


0% 


140 


0.771218(27) 


1.139(29) 


-2.39(32) 


5.0 


11 


93 


1% 


150 


0.771216(27) 


1.143(30) 


-2.44(34) 


4.8 


10 


90 


1% 


180 


0.771218(31) 


1.137(42) 


-2.36(54) 


4.8 


9 


85 


0% 


200 


0.771206(33) 


1.165(52) 


-2.79(70) 


3.9 


8 


86 


6% 


250 


0.771204(37) 


1.170(65) 


-2.86(95) 


3.9 


7 


79 


3% 


300 


0.771194(42) 


1.199(87) 


-3.38(1.41) 


3.6 


6 


72 


6% 


400 


0.771177(53) 


1.257(136) 


-4.50(2.46) 


3.3 


5 


64 


9% 


500 


0.771187(62) 


1.218(186) 


-3.68(3.65) 


3.2 


4 


51 


9% 


700 


0.771244(85) 


0.936(340) 


3.16(7.80) 


2.3 


3 


52 


2% 


1000 


0.771118 142) 


1.680 755 


-16.96(19.82) 


1.0 


2 


59 


7% 



Table 16: Fit (R 2 e ) = N 2v '(ae + beN' 1 + c e iV~ 3/2 ). DF is the number of degrees of freedom 
and CL is the confidence level of the fit. 



jVmin Og bg Cg X 2 DF CL 



60 


0.108209(3) 


0.107(2) 


0.08(2) 


17.8 


15 


27 


5% 


80 


0.108207(4) 


0.109(3) 


0.06(3) 


16.3 


14 


29 


3% 


100 


0.108204(4) 


0.113(4) 


0.02(4) 


14.3 


13 


35 


6% 


120 


0.108204(5) 


0.113(5) 


0.02(5) 


14.3 


12 


28 


5% 


140 


0.108206(5) 


0.110(6) 


0.05(6) 


13.5 


11 


25 


9% 


150 


0.108206(5) 


0.110(6) 


0.06(8) 


13.5 


10 


19 


6% 


180 


0.108208(6) 


0.106(8) 


0.11(11) 


13.0 


9 


16 


4% 


200 


0.108207(7) 


0.108(10) 


0.08(14) 


12.9 


8 


11 


7% 


250 


0.108205(8) 


0.113(13) 


0.00(19) 


12.5 


7 


8 


5% 


300 


0.108208(9) 


0.106(17) 


0.12(25) 


12.0 


6 


6 


1% 


400 


0.108203(11) 


0.120(26) 


-0.14(45) 


11.5 


5 


4 


2% 


500 


0.108201(13) 


0.128(37) 


-0.30(72) 


11.5 


4 


2 


2% 


700 


0.108212(17) 


0.076(64) 


0.95(1.46) 


10.5 


3 


1 


5% 


1000 


0.108185(30) 


0.227(156) 


-3.03(4.03) 


9.4 


2 





9% 



Table 17: Fit (R 2 g ) = N^iag + bgN" 1 + c g N^/ 2 ). DF is the number of degrees of freedom 
and CL is the confidence level of the fit. 
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N ■ 

1 * mm 




h 

u m 


Cm 


X 


171 


V^lj 


oU 




fl A A 7(1 A \ 
U.44 <\i-'±) 




1 Q R 
lo.O 


1 1 

LZ 


o no/ 
y.y /o 


oU 


n QQQ.ni q^oq^i 
U.ooyu loyZiO ) 


A 7Qf 1 8"\ 
U.4 ( y \lo) 


1 Q7f 1 
— l.o I yio ) 


1 1 Q 

1 l.o 


1 1 


/II 7°^ 
41. / /o 


1 nn 

1UU 








11. u 


1 n 

1U 


Q.^ 7% 


120 


0.339019(27) 


0.469(28) 


-1.23(34) 


9.8 


9 


36.4% 


140 


0.339038(29) 


0.436(33) 


-0.76(42) 


6.3 


8 


61.7% 


150 


0.339044(30) 


0.426(37) 


-0.61(48) 


5.9 


7 


55.6% 


180 


0.339041(32) 


0.431(41) 


-0.69(55) 


5.8 


6 


44.8% 


200 


0.339028(34) 


0.457(49) 


-1.09(69) 


4.9 


5 


43.1% 


250 


0.339010(41) 


0.496(68) 


-1.69(1.01) 


4.2 


4 


38.0% 


300 


0.339005(43) 


0.508(75) 


-1.88(1.13) 


4.1 


3 


25.5% 


400 


0.338976(62) 


0.583(140) 


-3.21(2.37) 


3.6 


2 


16.1% 


500 


0.338964(74) 


0.623(194) 


-4.04(3.62) 


3.6 


1 


5.9% 



Table 18: Fit (R 2 m ) = N 2u (a m + ^iV" 1 + c m N~ 3 / 2 ). DF is the number of degrees of 
freedom and CL is the confidence level of the fit. 



-^min 


a A 


b A 


ca 


2 

X 


DF 


CL 


60 


0.140305(6) 


-0.062(3) 


0.47(2) 


8.4 


15 


90.7% 


80 


0.140306(6) 


-0.064(5) 


0.49(4) 


8.1 


14 


88.2% 


100 


0.140303(7) 


-0.060(6) 


0.45(6) 


7.1 


13 


89.6% 


120 


0.140306(8) 


-0.064(8) 


0.49(9) 


6.6 


12 


88.2% 


140 


0.140306(8) 


-0.064(9) 


0.50(10) 


6.6 


11 


83.0% 


150 


0.140307(9) 


-0.065(10) 


0.52(12) 


6.5 


10 


77.0% 


180 


0.140309(10) 


-0.071(13) 


0.59(17) 


6.2 


9 


72.5% 


200 


0.140310(11) 


-0.072(16) 


0.62(22) 


6.1 


8 


63.4% 


250 


0.140307(12) 


-0.065(21) 


0.51(31) 


5.9 


7 


55.7% 


300 


0.140312(14) 


-0.078(27) 


0.72(42) 


5.3 


6 


50.3% 


400 


0.140310(18) 


-0.072(42) 


0.61(74) 


5.3 


5 


38.2% 


500 


0.140304(21) 


-0.049(59) 


0.12(1.15) 


5.0 


4 


28.9% 


700 


0.140306(27) 


-0.060(104) 


0.39(2.38) 


5.0 


3 


17.4% 


1000 


0.140297(47) 


-0.007(244) 


-1.03(6.34) 


4.9 


2 


8.6% 



Table 19: Fit A N = a a + 1 + caN 3//2 . DF is the number of degrees of freedom 
and CL is the confidence level of the fit. 
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N ■ 

1 * mm 


a B 


b B 


c B 


x 2 


DF 


CL 


60 


0.439619(31) 


-0.062(21) 


-0.02(17) 


12.8 


12 


38. 


1% 


80 


0.439587(34) 


-0.024(27) 


-0.45(27) 


8.3 


11 


68. 


3% 


100 


0.439578(38) 


-0.012(34) 


-0.61(37) 


8.0 


10 


63. 


2% 


120 


0.439596(41) 


-0.040(42) 


-0.24(49) 


6.6 


9 


67. 


5% 


140 


0.439620(43) 


-0.082(49) 


0.34(61) 


4.0 


8 


85. 


6% 


150 


0.439629(45) 


-0.098(55) 


0.59(70) 


3.5 


7 


83. 


4% 


180 


0.439626(48) 


-0.091(61) 


0.48(80) 


3.4 


6 


75. 


2% 


200 


0.439614(51) 


-0.068(73) 


0.12(1.01) 


3.1 


5 


68. 


5% 


250 


0.439590(61) 


-0.014(102) 


-0.72(1.52) 


2.5 


4 


63. 


6% 


300 


0.439583(65) 


0.004(118) 


-1.02(1.79) 


2.4 


3 


48. 


5% 


400 


0.439558(89) 


0.073(207) 


-2.27(3.55) 


2.3 


2 


31. 


9% 


500 


0.439528(106) 


0.174(284) 


-4.35(5.35) 


2.0 


1 


15. 


6% 



Table 20: Fit B N = a B + bsN' 1 + c B N' 3 ^ 2 . DF is the number of degrees of freedom 
and CL is the confidence level of the fit. 



-^min 


CLf 


b F 


CF 


x 2 


DF 


CL 


60 


0.000043(59) 


-0.039(39) 


1.30(33) 


11.5 


12 


49. 


0% 


80 


0.000105(65) 


-0.112(52) 


2.14(51) 


6.7 


11 


81. 


9% 


100 


0.000117(71) 


-0.129(65) 


2.36(71) 


6.5 


10 


76. 


8% 


120 


0.000087(77) 


-0.082(79) 


1.72(94) 


5.4 


9 


79. 


4% 


140 


0.000041(82) 


-0.004(93) 


0.64(1.17) 


3.0 


8 


93. 


5% 


150 


0.000026(86) 


0.023(104) 


0.24(1.34) 


2.6 


7 


91. 


7% 


180 


0.000037(90) 


0.002(116) 


0.54(1.54) 


2.5 


6 


87. 


.2% 


200 


0.000064(98) 


-0.052(139) 


1.38(1.94) 


2.0 


5 


85. 


3% 


250 


0.000110(116) 


-0.153(195) 


2.95(2.89) 


1.4 


4 


83. 


9% 


300 


0.000133(124) 


-0.205(220) 


3.80(3.32) 


1.2 


3 


76. 


2% 


400 


0.000167(174) 


-0.298(398) 


5.45(6.75) 


1.1 


2 


58. 


1% 


500 


0.000189(206) 


-0.372(545) 


6.98(10.23) 


1.0 


1 


30. 


7% 



Table 21: Fit F N = a F + bpN- 1 + cpN~ 3 / 2 . DF is the number of de erees of freedom 
and CL is the confidence level of the fit. 
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N ■ 

1 v min 


/ 


x 2 


DF 


CL 


60 


0.75(5) 


12.7 


14 


55.4% 


80 


0.79(7) 


12.2 


13 


51.3% 


100 


0.73(8) 


10.7 


12 


55.4% 


120 


0.65(9) 


6.7 


11 


82.0% 


140 


0.58(10) 


4.1 


10 


94.5% 


150 


0.57(11) 


4.0 


9 


91.2% 


180 


0.60(11) 


3.5 


8 


90.3% 


200 


0.64(12) 


2.7 


7 


91.5% 


250 


0.68(14) 


2.4 


6 


88.1% 


300 


0.69(15) 


2.3 


5 


80.2% 


400 


0.64(25) 


2.3 


4 


68.6% 


500 


0.69(51) 


2.3 


3 


52.0% 


1000 


1.18(1.17) 


2.0 


2 


36.1% 



Table 22: Fit F N (Rl) N = f. DF is the number of degrees of freedom and CL is the 
confidence level of the fit. 
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